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the  industry  and  a discussion  on  limitations  of  this 
database. 

The  harvest  patterns  of  PNG  crocodiles  are  discussed 
based  on  analysis  of  the  preceding.  No  trends  in  numerical 
harvest  were  found  for  either  species,  although  some  trends 
were  detected  in  the  frequency  of  harvest  composition  for  £. 
novaeauineae . Harvest  compositions  strongly  differed  among 
populations  for  most  within  year  intra-specific  comparisons 
for  each  species.  Despite  being  better  protected  by 
existing  legislation  with  respect  to  adult  size  classes,  C. 
corosus  populations  experienced  far  higher  harvest  mortality 

Models  are  presented  for  the  estimation  of  nesting 
female  crocodilian  size  from  clutch  characteristics.  The 


(pulse  model)  and  C.  novaeauineae  (flow  model)  as  well  as 
their  harvest  implications  are  discussed.  Reproductively 
active  female  alligators  from  widely  disjunct  harvested 

smaller  females  from  nesting.  Most  of  the  nesting  female 
crocodiles  were  smaller  than  the  size  protected  by  law  and 
were  thus  exposed  to  harvest,  which  presumably  led  to  a 
higher  turnover  rate  of  the  breeding  segment  and  a reduced 
production  of  hatchlings. 


Ontogenetic  variation  in  relative  growth  of  the  skull 
is  quantified  and  described  for  £.  novaeouineae.  An 
assessment  of  its  utility  for  the  prediction  of  the  sex  and 
size  of  individuals  is  made.  Allometric  coefficients  for 
estimation  of  size  are  given  and  models  for  gender 
prediction  through  parametric  discriminant  analysis  or 
nonparametric  binary  tree  classification  methods  are 
provided. 


CHAPTER  1 
PREFACE 


The  focus  of  this  dissertation  is  on  the  demographics 
of  selected  tropical  and  temperate  crocodilian  populations 
based  upon  inferences  drawn  from  several  sources.  These 
include  (1)  crocodile  skin  industry  trade  statistics  and 
associated  harvest  patterns,  (2)  the  comparative  gross 
morphometry  and  reproductive  mode  of  nesting  female 
crocodilians,  and  (3)  a quantitative  analysis  of  relative 
growth  of  the  crocodilian  skull.  Such  rather  diverse  topics 
were  selected  to  address  separate,  but  little  known,  aspects 
of  crocodilian  population  biology.  In  some  respects  they 
pose  similar  questions  (i.e.  - How  to  interpret  size  and  sex 
structure  of  populations  from  unknown  sex  hides  or  skulls, 
or  from  clutch  data  where  the  parenting  female  remains 
unseen?),  but  researched  from  different  methodologies.  As 
such  this  will  help  assess  how  well  each  section  complements 
or  refutes  the  findings  of  each  other  and  lends  toward  a 
holistic  approach  for  the  conservation  and  management  of 
these  resources.  These  studies  stem  from  10  months  of 
fieldwork  between  1980-1987  on  three  crocodilian  species:  21 
) working  with  Hew  Guinea 


d saltwater  (£. 
crocodiles,  plus  nine  months  divided  between  Louisiana  a 
Florida  researching  American  alligators  (Alligator 


In  Chapter  2,  I present  a summary  of  the  crocodile  skin 
industry  trade  statistics  for  PNG  from  1969-1980  and  review 
some  difficulties  associated  with  its  use.  These  years 
cover  the  period  from  the  inception  of  legislated  management 
of  the  crocodilian  resource  through  various  changes  toward 
the  development  of  the  industry  as  it  stood  for  what  is  the 

(roughly  80%  of  the  harvest)  and  £.  porosus  were  harvested 
during  those  years  for  which  the  annual  species  and  size 
compositions  of  harvests  from  11  geographic  areas  are 
presented.  Data  prior  to  1969  are  sparse  and  difficult  to 
interpret,  but  they  are  also  included  and  indicate  that  the 
total  harvests  were  substantially  greater  during 
pre-management  years.  Host  of  these  data  were  generated 
from  obscure  sources  contained  in  dusty  files  that  were  not 
well  safeguarded  while  continually  shifted  from  one 
government  department  to  another.  By  means  of  this  chapter, 
a permanently  accessible  reference  base  is  created  for 
management  and  conservation  authorities. 

The  synopses  found  in  Chapter  2 serve  as  the  database 
for  analyses  performed  in  Chapter  3,  which  deal  with  the 
observed  harvest  patterns  and  trends  of  those  data.  The 


data  from  Chapter  2 provide  an  indirect  measure  of  size,  but 
contain  no  information  on  the  sex  or  age  of  individuals.  In 
Chapter  3,  I tackle  the  thorny  problems  of  retrospectively 
assigning  hide  data  to  biologically  meaningful  size,  sex, 
and  age  classes.  These  transformations  are  then  used  to 
derive  harvest  mortality  rates  on  populations  of  each 
species  to  assess  the  impacts  of  harvesting  on  and  among 
them.  Limitations  in  the  use  of  such  data  are  discussed  and 


areas  for  further  research  are  suggested. 

In  Chapter  4,  I investigate  various  relationships  among 
gross  morphometries  of  nesting  crocodilians  and  how  those 
values  may  be  indicative  of  observed  reproductive  output 
(clutch  mass  and  clutch  size) . I also  examine  how  known 
clutch  parameters  may  be  useful  in  predicting  the  sizes  of 
unseen  parenting  female  crocodilians  and  thus  used  to 
profile  the  nesting  segments  of  breeding  populations. 

Within  this  framework,  I further  explore  the  comparative 
reproductive  modes  of  selected  tropical  and  temperate 


The  subject  of  Chapter  5 is  one  that  has  hitherto 
received  little  attention — a quantitative  analysis  of 
relative  growth  patterns  of  the  crocodilian  skull.  General 
ontogenetic  trends  of  crocodilian  cranial  growth  have  been 
outlined  by  Mook  (1921b),  Kalin  (193S,  1955),  and  Brongersma 


(1941).  Subsequently,  Iordansky  (1973)  developed  cranial 


indices  useful  for  interspecific  comparisons  of  relative 
growth  among  the  order  Crocodylia.  All  these  studies  were 
limited  by  the  lack  of  quantitative  assessment  of  growth 
with  respect  to  the  size,  sex,  and  age  of  individuals.  A 
major  advancement  in  this  field  occurred  in  the  mid-1970's 
when  Dodson  (1975)  explored  the  functional  and  ecological 
significance  of  relative  skull  growth  within  known  size 
American  alligators  (A . mississippiensisl . Here  too, 
however,  the  relationships  with  the  gender  and  age  of 
individuals  were  not  examined.  In  Chapter  5,  I extend 
previous  investigations  on  this  subject  area  by  including 
large  data  subsets  from  known  gender  animals  from  a species 
(£.  novaeguineae)  for  which  age  classes  are  known  or 
assigned.  Analytical  models  are  thereby  developed  to 
portray  ontogenetic  changes  in  relative  skull  growth  of  this 
crocodile  from  a multidimensional  character  framework.  That 
knowledge  is  used  to  assess  the  utility  of  functional 
patterns  of  relative  growth  as  a predictor  of  gender  and  age 

Demographic  parameters  of  crocodilian  populations  are 
often  poorly  understood  or  not  known  at  all.  The  ability  to 
monitor  and  interpret  demographic  trends  of  populations 
plays  a key  role  in  the  establishment  of  conservation 
stategies  and  tactics.  Knowledge  gained  from  these  studies 
will  better  enable  conservation  and  management  authorities 
to  determine  population  structural  parameters  of  crocodilian 


resources,  be  they  subject  to  harvest  or  accorded  full 
protection . 


CHAPTER  2 

CROCODILE  SKIN  INDUSTRY  TRADE  STATISTICS 


rural  commerce  throughout  lowland  areas  of  Papua  New  Guinea 
(PNG;  Fig.  2-1)  for  several  decades.  This  entails  the 


individuals  engaged  in  the  buying  or  exporting  of 
crocodile  products  from  PNG  must  be  licensed  by  the  Hildlife 


under  terms  of  the  1966  Crocodile  Trade  (Protection)  Act. 


repealed  in  the  Sepik  and  Madang  provinces  from  1972  to 
early  1975,  but  was  reinstated  thereafter.  Under  terms  of 


authorisation  from  the  Wildlife  Division  prior  to  the  export 
of  any  consignment  of  live  crocodiles  or  crocodile  skins. 
Additionally,  all  shipments  of  crocodilian  products  were 
made  subject  to  inspection  by  enforcement  officers. 

Following  the  independence  of  PNG  in  September  1975, 
concern  was  expressed  in  restructuring  the  crocodile  skin 
industry  whereby  national  interests  would  be  promoted  over 
those  of  expatriates  who  had  historically  been  the  main 
beneficiaries  of  the  trade.  Toward  this  end,  the  Food  and 
Agricultural  Organization  of  the  United  Nations  (FAO)  joined 
forces  with  the  PNG  government  in  1977  to  vitalize 
assistance  for  management  of  the  industry.  A major  shift  in 
the  focus  of  the  industry  occurred  in  1979  when  a national 
buying  scheme  was  established  to  encourage  the  sale  of  small 
live  crocodiles  to  large  scale  commercial  enterprises  for 
rearing  to  counter  the  economic  waste  in  the  marketing  of 
small  size  skins.  Subsequent  legislation  led  to  the 
introduction  of  a minimum  size  skin  limit  of  7”  (18  cm) 
belly  width  in  February  1961.  More  detailed  reviews  of  the 
history  of  crocodile  exploitation  in  PNG  and  government 
policies  toward  this  resource  may  be  found  elsewhere 
(Downes,  1971,  1973;  Bolton,  1980;  Bolton  and  Laufa,  1982; 
Hollands,  1987) . 


In  1980  I joined  the  FAO  monitoring  effort  on 
biological  investigations  of  PNG  crocodiles  and  compiled 
harvest  data  on  these  species  during  1980-1981.  My  initial 


review  of  Che  data  base  indicated  major  flaws  in  the 
collection  and  tabulation  process  since  forms  submitted  by 
buyers  and  exporters  were  ill-suited  for  the  compilation  of 
biologically  meaningful  data  and  hence  were  left  to 
accumulate.  My  objectives  in  presenting  these  data  are 
twofold:  (1)  to  provide  a permanently  accessible  reference 

and  biologists  involved  in  the  crocodilian  skin  industry  and 
(2)  to  illuminate  the  inherent  difficulties  encountered  in 
compiling  these  data  and  explain  how  they  were  circumvented. 
Other  factors  not  addressed  herein,  such  as  the  conversion 
of  commercial  belly  width  measurements  to  more  biologically 
meaningful  units,  size  class  sex  ratios,  and  the  assignment 

an  analysis  of  the  harvest  patterns  exhibited  by  these  data 
(CHAPTER  3)  . 


Methods 


Most  data  presented  herein  were  tabulated  from  the 
export  application  records  on  file  with  the  PNG  Wildlife 
Division  headquarters  in  Waigani,  Port  Moresby.  The  few 
remaining  data  are  from  published  accounts,  correspondence 
records  on  file  at  the  Wildlife  Division,  or  were  obtained 
directly  from  trade  sources.  Raw  data  were  compiled  from 
every  available  approved  application  for  each  licensed 


recorded  were:  exporter,  date  of  application,  origin  of 
purchases,  species,  skin  type  (belly  skin  or  horaback) , skin 
size,  and  number  of  skins  per  size  class.  Additionally, 
several  thousand  small  live  crocodiles  were  purchased  by 
government  stations  during  1979  and  1980  and  resold  to  large 

obtained  directly  from  government  sources. 

The  inclusion  of  buyer  quarterly  returns  as  a data  base 
was  rejected  on  the  grounds  that 


format  therein  w. 


not  explicit  and  therefore,  by  and  largely  useless  for 

and  followup  on  those  forms  was  evident.  Even  had  those 
records  been  complete  and  up-to-date  it  would  have  been 
hopelessly  confounded  by  the  occurrence  of  one  buyer 
reselling  his  stock  to  another  buyer  (rather  than  exporting 
directly)  which  would  introduce  an  immeasurable  factor  of 
duplication. 

A number  of  difficulties  arose  i 
compilation  of  raw  data.  £ 
in  which  they  were  resolved,  are  outlined  b 


Exporter  Identification 

Exporters  who  were  licensed  for  trade  in  the  crocodile 
skin  industry  during  the  years  1969-1980  are 


A-l,  Appendix  A.  However,  not  all  licensed  individuals 
subsequently  exported  crocodile  skins.  Furthermore, 
examination  of  the  then  used  exporters  file  registry  (since 
revised)  revealed  a disturbing  gap  in  the  sequential  order 
of  licensing  that  left  20  license  numbers  unaccounted  for. 

I presumed  that  those  numbers  were  not  in  fact  ever  used 
except  for  being  tentatively  assigned  to  unknown  applicants 
who  later  withdrew  their  applications  or  had  them  refused  by 
the  Management  Authority.  Additionally,  some  exporter  files 
were  discovered  to  be  missing  (part  1 of  CE  [Crocodile 
Exporter]  2,  CE4,  CE6  [part  1],  CE7  [part  1],  CE9,  CE12, 

CE16  [part  1],  and  CE45  [part  1])  and  never  accounted  for. 

Of  these  missing  files,  three  are  known  to  contain  the 
records  for  major  hide  industry  exporters:  CE2,  CE4  (both 
Papua  Region  traders)  plus  a New  Guinea  Region  trader,  CE45. 
The  years  1969  to  1972  are  affected  by  their  loss.  Data  for 
1973  onward  are  believed  to  be  a complete  reflection  of  all 
legal  exports. 

Pate  of  Export 

One  of  three  alternatives  was  available  for  this 
variable,  each  of  which  would  have  yielded  different  results 
on  a calendar  year  basis.  These  were:  (A)  the  exporter's 
date  of  application,  (B)  the  date  of  application  approval  by 
the  wildlife  Division,  and  (C)  the  date  of  intended  export. 


the  exporter' : 


application  in  tabulating 


data  as  that  most  nearly  correlated  with  the  time  of  actual 
harvest  and  afforded  the  minimum  inherent  overlap  between 
the  close  and  start  of  successive  calendar  years.  The  date 
of  export  approval  by  Wildlife  Division  was  not  used  because 
the  time  lag  bias  between  harvest  and  export  was  extended  by 
one  to  several  weeks  in  all  cases,  principally  due  to  mail 
delay  in  receipt  and  sometimes  on  account  of  irregularities 
in  filing  (form  incomplete,  lack  of  signature,  etc.). 
Likewise,  the  date  of  intended  export  offered  the  poorest 
correlation  with  the  time  of  actual  harvest  and  created  the 
greatest  overlap  between  years.  Furthermore,  the  date  of 
intended  export  was  often  not  specified  beyond  'earliest 
vessel'  or  else  altered  after  the  application  had  been 
approved.  Earlier  unpublished  in-house  attempts  at  export 
analysis  were  not  consistent  in  this  approach  and/or  made  no 
review  of  the  data  base  and  consequently  were  apocryphal . 

As  harvesting  occurred  throughout  the  year,  no 
corrections  for  annual  overlap  in  harvest  records  were  made 
since  all  years  were  so  affected  and  because  the  true 
measure  of  overlap  was  impossible  to  determine  (i.e.,  an 
application  dated  14  January  contained  skins  that  were 
primarily,  if  not  exclusively,  harvested  during  the  previous 
year) . As  a result  of  this  immeasurable  bias  and  because 
such  tabulations  would  prove  needlessly  unwieldy,  data  were 
not  compiled  to  reflect  a monthly  analysis  of  exports.  A 
priori  knowledge  indicated  that  hunter  effort  focused  on 


periods  of  little  or  no  moonlight  and  was  most  concentrated 
during  the  dry  season  when  crocodiles  were  more  accessible 
due  to  dwindled  aquatic  habitat  (Hall,  1983) . Marked 
variations  occurred  in  wet/dry  seasonality  between  harvest 


One  of  the  obvious  shortcomings  in  reviewing  data  from 
export  applications  was  that  the  origin  of  purchases  could 
only  be  identified  with  inordinate  difficulty  or  else  not  at 
all.  At  best,  the  origin  could  be  attributed  to  a 
particular  province  and  then  the  region  (Papua  or  New 
Guinea) , but  not  to  a specific  lake/lagoon  or  river  system. 
When  an  individual  application  specified  more  than  one 
province,  and  all  were  located  within  the  same  region,  that 
particular  consignment  was  accordingly  tallied  under  the 
category  of  Mixed  Papuan  Provinces  or  Mixed  New  Guinean 
Provinces.  When  an  export  permit  indicated  combined  skin 
purchases  from  both  regions  (i.e.,  Western  and  Sepik 
Provinces)  no  segregation  as  to  origin  was  possible.  In 
such  an  instance  the  only  recourse  was  to  record  the 
shipment  as  'unidentified  origin. ■ Papua  Region  provinces 
for  which  1969-1980  exports  were  identifiable  are:  Western, 
Gulf,  and  Central.  Those  so  known  from  the  New  Guinea 
Region  include:  East  and  West  Sepik,  Madang,  Morobe, 
Northern,  East  and  West  New  Britain,  and  the  North  Solomons 


(Fig.  2-1) . Of  necessity,  the  two  Sepik  provinces  and  the 

origin;  in  the  first  instance  this  was  because  the  provinces 
were  not  separate  political  entities  until  1975  and  in  both 
cases  this  was  due  to  inconsistent  identification  on  the 
export  permits.  However,  most  of  the  data  (probably  > 90%) 
tabled  as  Sepik  origin  stemmed  from  the  East  Sepik  Province. 


The  skins  of  the  two  crocodile  species  exported  from 

distinguished  from  each  other  by  a trained  observer.  New 
Guinea  crocodiles  possess  about  25  transverse  ventral  rows 
between  the  collar  and  anterior  vent  (see  Fig.  2-2),  while 
saltwater  crocodiles  have  30  or  more  rows.  Hence,  traders 


the  latter  species  as  'small  scale'.  Thus,  the  only 
consideration  here  was  the  implicit  assumption  that  an 
exporter  properly  listed  a consignment  under  the  appropriate 
species  and  skin  type  columns  when  filing  an  application. 

applications  as  to  the  species  identification  of  excessively 
large  individual  animals  (see  Tables  2-4  and  2-11) . 

Skin  Type  and  Size 


Skin  type  was  classified  either  as  belly  skin  (dorsal 
incision)  or  hornback  (ventral  incision) . The  methods  of 


measurement  for  each  type  are  legally  defined  in  PNG  (see 
Parker,  1981)  and  are  illustrated  in  Figures  2-2  and  2-3. 

It  should  be  noted  that  hornbacks  prior  to  February  1981  (as 
reported  herein)  were  formerly  measured  from  the  distal 
margin  of  the  abdominal  cut  to  the  proximal  edge  of  the 
third  row  of  dorsal  scutes  posterior  to  the  nuchals  on  each 
half  of  the  skin  and  then  summed  to  arrive  at  the  actual 
corresponding  belly  skin  size. 

Skins  and/or  live  crocodiles  were  categorized  into  size 
classes  based  on  one  inch  commercial  belly  width  increments 
as  specified  on  export  application  forms.  In  general,  this 
included  size  classes  ranging  from  4 to  20  inches  belly 

restriction  on  skin  size  from  Sepik  and  Madang  Provinces. 


Tallying  the  numbers  of  skins  per 
normally  a very  straightforward  affair 
two  exceptions. 


the  following 


First,  live  hatchlings  of  less  than  4 inch  commercial 
belly  width  were  purchased  during  1979  and  1980.  These  were 
included  within  the  4 inch  size  class  as  segregation  into 
actual  units  was  not  possible  or  practical  since  the  year 
class  was  identical. 

More  importantly,  between  1972-1979  one  or  more 
exporters  incorrectly  filed 


export  applications  by 


Figure  2-3.  Illustration  of  crocodile  hornback 
(modified  from  Parker,  1981)  - ' 

in  Papua  New  Guinea  


~je  distal  margin 


ww  legal  definition 
include  only  the  width  -* 
were  then  restricted  t< 


. - - - -w  hornback  skins 

sured  from  the  ventral  incision 
third  row  of  keeled  dorsal 
ed  so  as  to  exclude  the  dorsal 
Id  not  exceed  20"  (51  cm) . In 
hornbacks  was  changed  to 
31  tne  third  dorsal  row.  site  limits 
skins  between  2.5"  (6  cm)  and  6.5  " 


from  certain  size  classes 


combining  the  number  of  skins 
instead  of  reflecting  the  true  size  compositions.  The 
common  practice  of  those  buyers  was  to  purchase  skins  of 
4/5,  6/7,  and  8/9  inches  on  a same  value  per  piece  basis. 
These  skins  were  accordingly  lumped  on  export  applications 
which  effectively  reduced  the  number  of  size  classes  within 
that  range  from  six  to  three.  To  complicate  matters 
further,  a uniform  approach  in  recording  skins  was  not 
followed  by  these  exporters.  When  reporting  these  data  they 
sometimes  listed  them  under  the  lower  size  class  (e.g.,  6"), 

used.  This  fact  was  evidently  overlooked  or  disregarded  by 
earlier  workers  when  calculating  size  distributions  of  the 
harvest  and  average  skin  size.  I circumvented  this 
difficulty  by  segregating  all  lumped  data  found  on  export 
permits  and  adjusted  them  to  reflect  the  observed  summed 
ratios  for  those  size  classes  from  all  exporters  for  that 


The  trade  statistics  of  the  crocodile  skin  industry  in 
PNG  for  the  period  1969-1980  are  presented  in  a series  of 
tables  as  follows.  Composite  national  summaries  indicating 
the  political  regions  and  provinces  of  origin  for  each 
species  are  shown  in  Table  2-1  (£.  novaeguineae)  and  Table 
2-2  (£.  norosus) . Regional  harvest  summaries  reflecting 
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size  compositions  are  presented  in  Tables  2-3  - 2-4  for  £. 
novaeauineae  and  Tables  2-5  - 2-6  for  £.  porosus . Tables 
2-3  - 2-6  exclude  animals  of  unidentified  origin  listed  in 
Tables  2-1  - 2-2.  Annual  national  harvest  summaries 
indicating  size  compositions  by  provincial  origin  appear  in 
Tables  2-7  - 2-18  (£.  novaeguineae)  and  Tables  2-19  - 2-30 
(£.  porosus 1 . Several  thousand  small  live  crocodiles  were 
harvested  in  1979-1980  as  rearing  stock  for  large  scale 
commercial  ranches.  Those  animals  were  included  in  the 
previous  tables.  The  size  compositions  of  those  individuals 
are  presented  separately  in  Tables  2-31  - 2-34.  Data  for 

I was  unable  to  account  for.  Nevertheless,  total  harvest 
figures  (but  not  species  or  size  compositions)  were 
available  for  those  years  (Table  2-35) . 

Few  data  were  available  prior  to  legislative  attempts 
to  regulate  the  crocodile  skin  industry.  They  permitted 
even  less  interpretation  than  the  foregoing  and  are  shown  in 
Tables  2-36  and  2-37. 


Industry  Background 

Hollands  (1987)  identified  four  areas  of  concern  that 
historically  affected  the  crocodile  skin  industry  in  PNG: 
over-exploitation,  biologically  inefficient  harvesting, 
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for  national  harvest  of  Crocodvlus  novaeouineae.  Papua  New 
Guinea,  1970.  Data  incomplete  due  to  missing  files. 
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Table  2-19.  Size  compositions  (in  commercial  belly  width  inches) 
for  national  harvest  of  Crocodvlus  porosus.  Papua  New  Guinea, 
1969.  Mean  size  based  on  cohorts  4-20.  Data  incomplete  due  to 
missing  files. 


Western  Papua  Morobe  New  Guinea  Totals 

Prov.  Region  Prov.  Region 


Table  2-20.  Size  compositions  (in  commercial  belly  width  inches) 
for  national  harvest  of  Crocodvlus  porosus.  Papua  New  Guinea, 
1970.  Data  incomplete  due  to  missing  files. 


Western  Papua  Morobe  New  Guinea  Totals 

Prov.  Region  Prov.  Region 


Percent 


Table  2-21.  Size  composition  (in  commercial  belly  width  inches} 
for  national  harvest  of  Crocodvlus  porosus . Papua  New  Guinea, 
1971.  Data  incomplete  due  to  missing  files. 
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Table  2-31.  Live  crocodvlus  novaequlneae  purchases  included  in 


Prov.  Region 


Table  2-33.  Live  Crocodvlus  porosus  purchases  Included  in  1979 
national  harvest  summary  (Table  2-29) . Size  expressed  in 
commercial  belly  width  inches. 


Size  Western 


Table  2-34.  Live  Crocodvlus  porosus  purchases  included  in  1980 
commercial  belly  width  inches. 


Size  Western  Gulf  Central  Papua  Sepik  Hew  Guinea  Totals 
Prov.  Prov.  Prov.  Region  Prov.  Region 
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economics,  and  poor  international  marketing.  A brief 

Prior  to  the  imposition  of  legislative  controls,  the 
preferred  trade  in  PNG  was  for  large  adult  crocodiles 
(especially  £.  porosusi  which  rapidly  diminished  the  wild 
breeding  stock.  Consequently,  hunters  increasingly  relied 
on  the  cropping  of  small  animals  with  little  economic  value. 
This  resulted  in  dramatically  larger  harvests  during  the 
1960 's  and  altered  population  compositions  to  the  extent 
that  conservationists  feared  the  harvest  could  no  longer  be 
sustained  and  that  PNG  crocodiles  were  imperiled  with 
extinction  (Neill,  1971).  Fortunately,  those  views  proved 
to  be  rather  alarmist.  Nevertheless,  they  did  help  inspire 
early  management  actions. 

Concurrently,  villager  hunters  derived  marginal  returns 
on  the  true  value  of  their  resource  from  expatriate  traders. 
Additionally,  national  revenues  were  suppressed  since  the 
bulk  of  exports  were  marketed  through  Singapore  middlemen, 
rather  than  through  direct  contact  with  European  and 
Japanese  tanneries.  A further  combination  of  poor  skinning 
and  preservation  techniques,  plus  the  inclusion  of  less 
valuable  homback  hides,  helped  to  depress  the  potential 
economic  return.  Trade  in  hornback  hides  (principally  from 
the  Sepik  provinces)  was  very  popular  with  the  Japanese 
market  in  the  mid-1970's  (it  comprised  nearly  half  of  the 
1973  total  exports) , but  diminished  to  almost  nil  by  the  end 


of  the  decade.  A colorful  perspective  on  the  life  and  tines 
of  a crocodile  trader  during  this  period  is  provided  by 
Pinney  (1976) . 

During  1976  a marketing  service  for  crocodile  skins  was 
established  by  the  national  government  to  help  ensure  fair 
market  competition  with  private  sector  interests.  This 
action,  along  with  extension  activities  aimed  at  improving 
the  quality  of  hides  through  better  skinning  and 
preservation,  greatly  enhanced  revenue  return  among  village 

conservation  and  management  policies.  These  latter  aspects 
were  further  strengthened  through  a policy  shift  in  the  late 
1970's  to  replace  the  hunting  for  skins  of  very  small 
crocodiles  by  the  live  harvest  of  those  animals  for  rearing 
to  commercial  size.  Under  this  policy  hunters  received  far 
higher  prices  for  a live  animal  than  for  a comparable  size 
skin  and  thus  could  earn  a greater  income  from  harvesting 
fewer  crocodiles.  X refer  readers  to  Hollands  (1987)  for  a 
more  thorough  authoritative  treatment  on  this  subject. 

The  otherwise  excellent  review  paper  by  Hollands  (1987) 
on  the  management  of  Papua  New  Guinea  crocodiles  is  marred 
by  inaccurate  harvest  figures  (p.  76)  for  years  in  which 
data  could  be  substantiated,  rather  than  speculated.  Table 
2-38  indicates  those  data  along  with  the  comparative 


Hollands  (19 
£.  n-  £.  E-  T 


corrected  values  reported  herein.  In  some  instances  the 
differences  appear  to  be  minor,  but  in  one  instance  (the 
1974  Q.  novaeouineae  harvest)  there  is  an  egregious  doubling 
of  the  actual  value.  I assume  that  the  figures  reported  by 
him  for  the  1968-1969  to  1972-1973  period  represent  fiscal 
year  (July-June)  totals.  From  1973  on  he  used  calendar  year 
exports,  as  I have  throughout  this  study.  The  pre-1969 
harvest  figures  in  that  paper  are  admittedly  speculative  and 
their  degree  of  accuracy  is  impossible  to  ascertain, 
although  they  do  not  seem  unreasonable  to  me. 

Similarly,  the  harvest  figures  presented  by  Whitaker 
(1980:45-46)  and  Whitaker  and  Kemp  (1981:23  and  25)  are 
largely  inaccurate.  They  were  based  on  preliminary 
tabulations,  after  which  additional  records  were  available. 
The  correct  values  for  the  appendices  in  those  reports  may 
be  obtained  from  the  data  in  this  study. 

There  are  several  other  sources  of  PNG  crocodile 
harvests  not  reflected  in  this  study.  From  1969-1975  an 
unknown  number  (2,000-3,000?)  of  crocodiles  (mostly  £. 
novaeouineae)  of  unknown  population  structure  were  caught 
live  in  Western  Province  and  sold  to  government 
demonstration  farms.  From  1976-1980  a further  unknown 
number  (5,000-10,000)  of  primarily  juvenile  C.  novaeouineae 
were  caught  live  and  placed  on  village  ranches  in  Western 
Province.  This  also  occurred,  to  a lesser  extent,  in  Gulf 
Province  and,  perhaps  to  a 


Province.  During  that  same  period,  compensation  claims  of  5 
kina  (about  US$7)  were  paid  to  villagers  in  East  Sepik 
Province  for  oversize  crocodiles  (>  20"  belly  width) 
inadvertently  or  intentionally  killed.  The  total  number  and 
composition  of  animals  thus  'bountied'  could  not  be 
obtained,  but  I believe  it  to  be  between  500-1,500.  This 
harmful  practice  was  eventually  discontinued  in  the  early 


Although  data  for  the  years  1973-1980  are  complete, 
only  about  60%  of  the  1972  data  are  represented,  while  data 
for  1969-1971  were  largely  unaccounted  for.  Data  for 
earlier  years  were  less  accessible  still,  and 
correspondingly  more  difficult  to  quantify.  Records  were 
shifted  among  several  governmental  departments  over  the 
years  and  their  maintenance  suffered  as  a consequence. 

The  size  class  compositions  shown  herein  are  assumed  to 
be  representative  of  the  populations  at  large  with  the 
following  exceptions.  Bias  exists  in  the  number  of  animals 
> 20"  belly  width  harvested  during  1972-1975  since  illegal 
transhipment  of  oversize  animals  from  Western  Province  to 
the  Sepik  provinces  occurred  (pers.  comm,  from  pilots  and 
Lake  Murray  hunters) . Second,  the  largest  size  class 
thereafter  (20")  is  disproportions! ly  represented  through 


inclusion  of  oversize 


shrunken 


legal  requirements.  Third,  a marketing  bias  existed  against 
skins  < 8"  belly  width  and  they  are  thus  unrepresented.  For 
example,  three  buyers  handled  > 95%  of  Western  Province 
skins  during  1979:  Daru  Trading  (CE81),  Laki  Trading  (CE71) , 
and  Wildlife  Division  (CE65) . As  a matter  of  policy,  CE81 
did  not  purchase  skins  < 7",  CE  bought  live  crocodiles  but 
not  skins  < 8",  and  < 5%  of  CE71's  Western  Province  records 
(N  - 4,678)  were  < 6".  Lastly,  hunter  selection  against 
harvesting  small  crocodiles  is  presumed  to  have  occurred  in 

rearing/holding  facilities  did  not  exist. 

All  but  a negligible  portion  of  the  data  in  this  study 
(<0.1  %)  stem  from  the  harvest  of  wild  animals,  although  the 
production  of  crocodile  skins  from  commercial  ranches  now 
forms  a significant  part  of  the  annual  trade  in  PNG.  While 
approximately  80  exporters  were  known  to  be  licensed  during 
the  1969-1980  period  (Table  A-l,  Appendix  A) , the  bulk  of 
the  trade  in  any  given  year  lay  in  the  hands  of  only  a few 
licensees.  In  fact,  by  1980  a single  license  holder  was 
responsible  for  more  than  40%  of  the  total  exports. 

Several  important  events  occurred  in  response  to  the 
compilation  of  these  data.  The  difficulties  encountered 
provided  the  impetus  for  the  establishment  of  a secure 
computerized  data  base  for  these  and 


future  records  within 


the  Wildlife  Division.  They  also  drew  attention  to  the 
inadequacy  of  existing  documentation  of  the  trade. 
Accordingly,  standardized  purchase  dockets  were  designed  and 
subsequently  issued  free  to  each  export  license  holder. 

These  forms  better  addressed  the  biological  needs  of 
management,  yet  simultaneously  allowed  confidentiality  of 
financial  transactions  desired  by  licensees.  Furthermore, 
the  harvest  record  summaries  of  this  study  provided 
persuasive  quantitative  evidence  of  the  economic  wastage  in 
the  trade  of  very  small  crocodile  skins  and  contributed  to 
the  introduction  of  a lower  size  limit  in  19B1. 

Following  the  successful  shift  of  management  policy 
toward  the  purchase  of  small  live  crocodiles  for  rearing  and 
the  establishment  of  large  scale  commercial  ranches  during 
the  late  1970's  and  early  1980's,  FAO  terminated  its 
involvement  in  the  PNG  crocodile  project.  A Crocodile 
Advisory  Board,  consisting  of  various  government  and 
industry  representatives,  was  formed  to  help  fill  that  void. 

Within  little  more  than  a decade  the  crocodile  skin 
industry  in  PNG  was  transformed  from  a troubled,  unregulated 
industry  with  a suspect  future  to  one  of  a competently 
managed  resource  that  continues  to  assure  the  economic 
livelihood  of  many  rural  people  while  maintaining  healthy 
populations  of  both  crocodile  species. 


GUINEA, 


For  several  decades  the  New  Guinea  (Crocodvlus 


saltwater  (£.  porosus)  crocodiles  have 


susceptible  to  capture  (Hall,  1983) . 


commercial  belly  width  (see  METHODS)  was  imposed  on  skin 
trafficking  for  both  crocodile  species.  Political 
opposition  led  to  the  repeal  of  the  51  cm  clause  in  the 
Sepik  and  Madang  provinces  until  early  1975,  when  it  was 
reinstated,  since  that  time,  the  upper  size  limit 
restriction  has  been  uniform  throughout  PNG.  Additionally, 
under  the  1966  law  skin  shipments  were  made  subject  to 
inspection  by  enforcement  officers.  Further  still,  a 
requisite  of  licensing  stipulated  that  buyers  file 
individual  quarterly  returns  on  their  purchases  and  that 
exporters  first  secure  written  authorization  from  the 
Wildlife  Division  prior  to  export  of  any  consignment  of  live 
crocodiles  or  crocodile  skins.  Subsequent  legislation  led 
to  the  introduction  of  a minimum  size  skin  limit  of  18  cm, 
effective  as  of  15  February  1981  (Parker,  1981) . Other 
reviews  of  the  history  of  crocodile  exploitation  in  PNG  and 
government  policies  toward  this  resource  have  been  presented 
elsewhere  (Downes,  1971,  1973;  Bolton,  1980;  Bolton  and 
Laufa,  1982;  Hollands,  1987). 

As  part  of  an  FAO  (Food  and  Agricultural  Organization 
of  the  United  Nations)  monitoring  effort  on  biological 
investigations  of  PNG  crocodiles,  I compiled  harvest  data  on 
these  species  during  1980-1981  (CHAPTER  2) . My  initial 
review  of  the  data  base  indicated  major  flaws  in  the 
collection  and  tabulation  process.  Forms  submitted  by 
buyers  and  exporters  were  ill-suited  for  the  compilation  of 


biologically  meaningful  data  and  hence  were  left  to 
accumulate  without  being  analyzed  once  shipments  were 
inspected  and  approved.  Some  of  these  problems  have 
subsequently  been  overcome  (Hollands,  1984,  1987).  The 
present  data  set  is  an  interpretation  of  those  harvests 
prior  to  reforms  initiated  to  ensure  an  adequate  information 
retrieval  system,  and  it  highlights  some  of  the  difficulties 
associated  with  such  data. 

My  objectives  in  examining  these  data  were:  (1)  to 
determine  the  size,  composition,  and  origin  of  PNG  crocodile 
harvests  by  species;  (2)  to  detect  whether  patterns  and 
long-term  trends  occurred  within  and  among  populations  of 
each  species;  and  (3)  to  assess  the  biological  and  economic 
impacts  of  harvesting  on  the  crocodile  resource.  In 
essence,  this  research  sought  to  learn  if  existing 
legislation  and  policies  were  practicable  and  if  they 
provided  adequate  safeguards  to  ensure  sustainable  yields  to 
rural  people  who  relied  on  crocodiles  as  a basis  for  their 


Harvest  data  for  1969-1980  were  compiled  from  licensed 
dealer  export  permits  on  file  at  Wildlife  Division 
headquarters  in  Waigani,  National  Capital  District  (9'30'S, 
147'11'E).  Data  were  tabulated  on  a calendar-year  basis  by 
the  exporters  date  of  application  as  that  most  closely 


paralleled  the  time  of  harvest.  Variables  recorded  were: 
exporter,  date  of  application,  provincial  origin  of 
purchases,  species,  skin  type  (dorsal  incision  = belly  skin, 
ventral  incision  - hornback) , skin  size,  and  number  of  skins 
in  each  size  class.  Hide  counts  were  tabulated  into  size 
classes  as  recorded  by  commercial  belly  width  (BW,  which  is 
approximately  70%  of  the  chest  girth;  see  Parker,  1981) . 
These  original  data  were  grouped  in  one  inch  (2.54  cm)  BW 
increments  as  specified  by  existing  regulations,  and  ranged 
from  4-20  inches  (10  - 51  cm)  BW.  In  this  study,  a 
correction  of  4.4%  was  made  for  observed  shrinkage  in  hides 
(N  ■ 34,  SE  » 0.46)  to  derive  preskinning  belly  width  size 
classes,  which  were  then  transformed  to  equivalent  total 
length  (TL)  measurements.  Transformations  for  C. 
novaeouineae  follow  Montague  (1984b;  equations  10-12,  39-40; 

not  13-50  cm  - Montague,  in  litt.J).  The  conversion  of 
belly  width  to  total  length  measurements  for  PNG  c.  porosus 
follows  Whiteside  (1978;  y = 3.5  + 4.05  x).  Prior  to  1969, 
few  data  were  available  for  trade  in  crocodile  skins,  and 
they  permitted  only  limited  interpretation. 

In  addition  to  skins,  live  crocodiles  were  purchased  by 
Wildlife  Division  demonstration  crocodile  farms  for  resale 
as  stock  to  large-scale  commercial  crocodile  ranches.  A 
breakdown  of  live  crocodile  purchases  in  1979  and  1980  was 
obtained  from  government-station  purchase  dockets  issued  at 


Angoram  (4*01'S,  144‘01'E;  East 
(7‘01'S,  141'30'E;  Western  Province),  Balimo  (8*04'S, 

142 *59 'E;  Western  Province),  Green  River  (3*54'S,  141‘11'E; 
West  Sepik  Province),  Kikori  (7*27'S,  144*14'E;  Gulf 
Province),  Hoitaka  (9*30'S,  147*11'E;  Central  Province),  and 
Pagwi  (4‘01'S,  143‘0'E;  East  Sepik  Province). 

Sex  determination  of  hides  was  not  possible  since 
crocodilians,  except  some  adult  male  gharials  (Gavialis 
ganoeticus) , do  not  exhibit  secondary  sexual  characteristics 
other  than  size  dimorphism  (Rao,  1981) . Size  dimorphism  is 
not  apparent  in  PNG  crocodiles  under  2.7  m TL  and  3.5  m TL 
for  £.  novaeouineae  and  £.  porosus.  respectively.  Thus, 
except  for  a few  individuals  listed  in  Table  3-11,  size 
dimorphism  cannot  be  used  as  a basis  for  gender 
identification  of  the  data  in  this  paper.  Data  on  £. 
novaeouineae  sex  ratios  were  collected  by  cloacal 
examination  (Chabreck,  1963)  of  caught  wild  animals  placed 
in  rearing  pens  at  Angoram  in  East  Sepik  Province  (northern 
population)  and  in  two  subdistricts  of  Lake  Hurray  District 
in  Western  Province  (southern  population;  Fig.  3-2) . The 
latter  stock  were  sampled  at  Kaviananga  (7*36'S,  141’18'E) 
and  Komovai  (7*33'S,  141*15'E)  in  the  Middle  Fly  subdistrict 
and  at  Boikumava  (7*23'S,  141*23'E),  Egiza  (6*51'S, 

141*30'E) , Kapikam  (6*43'S,  141*45'E),  Magipopo  (6*50'S, 
141*27'E),  Tagum  (7*04'S,  141*34 'E) , Upovia  (6*40'S, 

141'27's) , and  Wamak  (Boimbulavu  [6*50'S,  141*26'E)) 


142°  itfE 


villages  in  the  Lake  Murray  subdistrict,  including 


government  crocodile  farm  at  Baboa.  Early  Baboa  data  were 
compiled  from  station  files.  Due  to  a paucity  of  data,  sex 
ratios  were  estimated  for  £.  novaequineae  between  177-214  cm 
TL  (Table  3-3)  and  for  all  £.  porosus  up  to  218  cm  TL.  The 
sex  ratios  of  these  species  were  arbitrarily  assumed  to  be 
similar,  and  mean  value  estimates  for  BW  size  class 


subgroups  (4-6",  7-10",  11-14" 
species  were  applied  to  the  latter  species  (Table  3-4) . 

This  was  essential  in  order  to  calculate  the  sex-based  age 
class  distributions  for  further  analysis  (Tables  3-6  and 
3-9)  following  the  assignment  of  sizes  to  age  classes  (Table 
3-5) . That  methodology  is  detailed  in  Appendices  B and  C. 

The  estimated  age-size  relationships  for  £. 
novaequineae  (Table  3-5)  are  based  on  my  interpretation  of 
unpublished  growth  data.  Growth  data  for  PNG  £.  porosus 
were  lacking  and  the  age-size  relationships  t 
(Table  3-5)  were  derived  from  a model  I 
porosus  (Webb  et  al.,  1978). 

Harvest  mortality  tables  (Tables  3-6  and  3-9) , 
patterned  after  dynamic  composite  life  tables  (Caughley, 

1977 t Seber,  1982),  were  constructed  from  pooled  harvest 
data  where  the  number  dying  per  age  class  (dx)  represents 
actual  summed  harvest.  Those  data  were  used  to 
back-calculate  the  number  alive  at  the  beginning  of  each  ag 
class  (lx)  and  to  estimate  overall  hunting  mortality  (qx) . 


r Australian  £. 


For  certain  analyses,  the 
e size  classes:  juveniles  ( 


were  subdivided  i 


(1.5 


class  structure  was  expressed  by  proportion 
he  immeasurable  interaction  of  factors  that 
d numerical  harvest.  Binomial  proportions  tests 
(Z-statistic;  Mendenhall,  1968)  were  used  to  test  the  null 
hypothesis  that  the  proportions  of  juveniles,  subadults,  and 

comparisons  (Tables  3-7  and  3-10) . In  this  instance  and 
elsewhere  (size  class  trend  analysis)  the  significance 
levels  (P  values)  were  adjusted  by  using  Bonferoni's 
inequalities  (Snedecor  and  Cochran,  1980)  because  the 
proportions  were  interdependent. 

I used  the  Mann-Kendall  test  (Bradley,  1968)  to  test 
for  the  existence  of  trends  (increase  or  decrease)  in 
numerical  harvests  and  the  Kendall  Tau  test  (Sokal  and 
Rohlf,  1973)  to  detect  trends  among  the  frequency  of  harvest 
size  classes.  Where  trends  were  detected,  I regressed  the 
size  class  proportions  against  the  time  of  harvest  (year 
number)  to  obtain  the  predicted  proportions  and  their 
respective  residuals.  The  number  of  positive  and  negative 
residual  values  were  then  tabled  against  successive  year 
numbers  of  harvest  and  analyzed  by  chi-square  contingency 


tests  (2x2,  2x3,  2X4)  to  evaluate  if  the  data  showed 
cyclic  periodicity. 


Total  harvest  figures  were  unobtainable  for  most  years. 
Files  for  the  years  1969-1971  were  incomplete  on  a 
nation-wide  basis  (either  lost  or  destroyed) . A similar 
situation  existed  for  1972  data  from  the  Sepik  Provinces 
only.  Complete  files  existed  for  the  remainder  of  the 
period,  but  from  1973  to  I960  large  portions  of  the  national 
harvest  were  not  identifiable  as  to  province  of  origin  for 
certain  years  (Tables  3-1  and  3-2) . 

Although  few  data  from  1969-1972  were  available  for 
analysis,  the  total  crocodile  harvest  for  that  period  was 
known,  though  it  was  not  identifiable  as  to  species,  origin, 
or  size.  Crocodile  exports  totaled  27,094  and  29,868  for 
1969  and  1970,  respectively  (Wildlife  Division  file  14-9-19, 
part  l),  while  exports  for  1971  and  1972  were  21,135  and 
45,185  (Post  Courier  News,  Port  Moresby,  PNG;  2 June  1973). 

Similarly,  some  data  were  available  for  years  prior  to 
1969,  but  the  species  and  size  distributions  were  not  known. 
A single  exporter  (Booth  & Co.)  purchased  a total  of  112  505 
PNG  crocodile  hides  during  1959-1968  (Wildlife  Division  file 
14-9-19,  part  1) . This  figure  provides  some  insight  into 
the  magnitude  of  unregulated  harvest  (pre-1969)  as  it  is 
equal  to  one-fourth  of  the  total  for  all  licensed  exports 
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from  1969-1980  (445,776).  The  provincial  origin  of  these 
pre-1969  hides  was  as  follows:  Milne  Bay  - 2,318:  Western  - 
51,488:  Westem/Gulf  - 26,005;  Gulf  - 2,319;  Central  - 2,823 
(Papua  Region  subtotal  = 84,953);  Sepik  - 27,243;  Manus  - 
309  (Mew  Guinea  Region  subtotal  - 27,552). 


c.  novaeguineae 

Export  file  records  for  PNG  harvests  of  £.  novaeguineae 
totaled  281,027  (annual  X = 31,225)  for  1972-1980,  the  years 


for  which  data  were  largely  complete.  Of  this  total,  50.9% 
(annual  x = 15,899)  was  from  the  Papua  (southern)  Region  and 
38.8%  (annual  x = 12,127)  from  the  New  Guinea  (northern) 
Region,  with  10.2%  of  unknown  origin  (Table  3-1).  Since 
Western  and  the  Sepik  Provinces  each  comprised  more  than  80% 
of  the  harvest  identifiable  as  to  origin  in  their  regions 
(Papua  and  New  Guinea,  respectively) , data  for  £. 
novaeguineae  from  those  provinces  were  considered  to  be 
broadly  representative  and  were  used  for  geographic 
comparisons  (Figs.  3-3  and  3-4). 

Excluding  hides  of  unknown  origin,  these  data  indicated 


a minimum  mean  annual  cro 
animals  with  an  average  s 
Based  on  discussion  with 
that  total  was  believed  t 
District  (Fig.  3-2) ; most 
Aramia-Bamu  River  basins 
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£.  novaeguineae  £.  porosus 


annual 


(Fig.  3-3) . Mean 
Western  Province  for  the  period  1956-1970  was  estimated  at 
7,000  (Whitaker,  1980).  Analysis  of  1966-1967  data 
(regarded  as  the  peak  of  historical  exploitation;  from 
Bustard,  1967)  revealed  a mean  skin  size  of  181  cm  TL,  (39.1 
cm  BW) , with  16.2%  of  the  sample  > 2.2  m TL.  The  average 
skin  size  for  that  interval  was  presumed  to  have  been  larger 
than  that  for  the  more  recent  period  due  to  a selective  bias 
by  hunters  and  buyers  against  < 0.8  and  > 2.2  m TL  size 
classes  during  1969-1980.  When  those  size  classes  are 


Data  for  the  Sepik  Provinces  showed  a larger  numerical 
harvest  than  for  Western  Province;  however,  the  size  at 
harvest  was  very  similar.  The  minimum  mean  annual  harvest 
for  1972-1980  was  8,903  with  an  average  skin  size  (excluding 
animals  > 2.2  m TL)  of  135  cm  TL  (21.0  cm  BW) . Mean  sample 


classes),  with  only  1.1%  > 2.2  m TL. 

Hunting  mortality  in  Western  Province  (20.2%)  was 
significantly  lower  (Z  = -34.9,  P < 0.001)  than  that  in  the 
Sepik  Provinces  (23.5%)  (Table  3-6). 

Harvest  data  were  complete  for  the  years  1973-1980  and 
no  trend  in  numerical  harvest  was  detected  for  total 
national  harvests  in  that  period  (Mann-Kendall  test,  P > 
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0.05).  Regional  numerical  harvests  could  not  be  so 
evaluated  due  to  the  large  numbers  of  animals  of 
unidentified  origin. 

However,  when  the  frequencies  of  harvest  size  class 
compositions  were  examined,  some  trends  were  detected  in 
Western  Province  during  the  1969-1980  interval.  Juveniles 
showed  an  increase  in  proportional  harvest  composition  while 
the  adult  size  class  experienced  a decrease  (Kendall  Tau 

Comparative  data  among  populations  existed  only  for  the 
years  1972-1980  (Fig.  3-5) . During  that  period  no  trends  in 
harvest  composition  were  detected  for  Western  Province 
harvests,  while  a decreasing  trend  in  the  subadult  size 
class  was  detected  in  harvests  from  the  Sepik  Provinces 
(Kendall  Tau  test,  P < 0.05).  Cyclic  periodicity  was  absent 
in  all  but  the  Western  Province  juvenile  size  class  for 
which  a 3 year  pattern  existed  (X2  = 6.00,  2 df,  P < 0.05). 
Significant  differences  in  size  class  proportions  were  found 
in  most  within-year  comparisons  among  populations  (Table 

The  proportion  of  male  to  female  juvenile  and  subadult 
C.  novaeouineae  (H  = 1,777)  in  Lake  Hurray  District,  Western 
Province  (Table  3-8)  was  the  same  among  sample  years  (Z  = 
0.48,  P > 0.05)  with  a 1:1  sex  ratio  for  1969  (Z  = 0.41,  P > 
0.05)  and  1980  (Z  - 0.45,  P > 0.05) . However,  both  the  Lake 


Murray  (male  biased)  and  Middle  Fly  (female  biased) 
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subdistrict  samples  differed. 


as  from  each  other  (Z  = 4.89,  P < 0.001).  Meanwhile,  a 
female  bias  was  observed  for  a small  sample  of  the  northern 
population  at  Angoram  on  the  Sepik  River  (Z  » -2.87,  P < 
0.005).  No  data  were  available  for  adult  £.  npyaegaineae - 
Although  the  composite  sex  ratio  for  Lake  Murray  £. 
novaeouineae  indicated  the  existence  of  an  equal  proportion 
of  males  to  females,  there  were  generally  significant 
departures  from  unity  within  size  classes  (Table  3-3) . 


£.  corosus 

Total  export  file  records  for  PNG  harvests  of  £. 
porosus  from  1972-1980  was  66,164  (annual  x - 7,352);  52.8% 
originated  from  the  New  Guinea  Region  (annual  x «=  3,883), 
2.5%  from  the  Island  Region  (annual  x - 435)  origin,  34.4% 
from  the  Papua  Region  (annual  x = 2,531),  and  10.1%  of 
unknown  locations  (annual  x = 1,117)  (Table  3-2).  The 
average  size  at  harvest  (excluding  animals  > 2.2  m TL)  was 
83  cm  TL  (19.6  cm  BW)  in  the  New  Guinea  Region,  113  cm  TL 
(27.0  cm  BH)  in  the  Island  Region,  and  81  cm  TL  (19.1  cm  BH) 
in  the  Papua  Region. 

Hunting  mortality  rates  in  the  Papua  Region  (59.7%)  and 
New  Guinea  Region  (56.6%)  were  similar,  but  statistically 
different  (Z  - 9.64,  P < 0.001).  However,  hunting  mortality 
was  significantly  lower  than 


Island  Region  (35.0%) 


either  the  Papua  Region  (Z  ■ 32.4,  P < 0.001)  or  New  Guinea 
Region  (Z  ■ 28.8,  P < 0.001)  (Table  3-9).  Hunting  mortality 
for  combined  £.  oorosus  populations  (56.8%)  was 
significantly  greater  than  that  for  C.  novaeouineae  (21.7%) 
(Z  = 241.5,  P < 0.001) . 

As  with  £.  novaeouineae.  no  trend  was  found  in  C. 
porosus  for  the  1973-1980  numerical  national  harvests 
(Mann-Kendall  test,  P > 0.05).  Unlike  the  case  for  the 
former  species,  no  trends  or  cyclic  patterns  were  detected 
in  C . porosus  size  class  compositions  (Kendall  Tau  test,  P > 

within-year  comparisons  among  populations  revealed 
significant  differences  in  harvest  size  class  proportions 

Few  data  (N  = 53)  were  available  for  juvenile  and 
subadult  £.  porosus.  but  they  indicated  a 0.6:1  male-female 
sex  ratio  (Z  - -3.03,  P < 0.005)  in  Lake  Hurray  District  of 
the  Papua  Region  (Table  3-8) . 

Harvesting  of  all  sizes  of  animals  of  both  species  was 
practiced  in  the  New  Guinea  Region  until  1975.  Partial  1972 
data  were  available  as  well  as  complete  data  for  1973-1975. 
Estimates  of  size/age  class  and  sex  ratios  were  not  derived 
for  crocodiles  > 2.2  m TL  to  obtain  estimated  harvest 
mortality  rates  for  such  individuals,  but  such  data  indicate 
that  minimum  totals  of  1,410  £.  ngyaequlnsae  (X  = 240  cm  TL) 
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harvested 


in  that  period  (Table  3-11) . 


Harvest  Trends 

Total  annual  harve 
reflection  of  hunting  a 


: of  PNG  crocodiles  is  largely  a 
I marketing  (domestic  and 


d wide  variations  existed  among 
level,  by  availability  of  salt. 


Figures  3-7  and  3-8 
ranges  in  precipitation 


international)  conditions 
and  within  populations  for  these 
strongly  influenced  by  water  le\ 
batteries,  electric  torches,  anc 
access  to  skin  buyers  (Hall,  19t 
illustrate  the  within  and  among 
patterns  that  affected  crocodile  harvests. 

The  observed  trends  within  and  among  these  crocodile 
populations  are  artifacts  of  selective  harvesting  to  satisfy 
market  demand  as  well  as  natural  dynamic  flux.  Such 
harvesting  entailed  selection  for  many  thousands  of  very 
small  (0.4  - 0.9  m TL)  juveniles  (primarily  in  the  New 
Guinea  Region  populations) , which  were  utilized  as  watch 
straps  in  the  Singapore  and  Japanese  markets.  Meanwhile, 
selection  against  those  sizes  of  animals  occurred  in  the 
Papua  Region.  Voluntary  conservation  measures  and  marketing 
logistics  in  the  latter  region  were  responsible  for  that 
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Figure  3-7.  Means  (horizontal  lines)  of  monthly 
precipitation  at  Bossett,  Papua  New  Guinea,  1965-1980  (no 
1968  data).  Rectangles,  standard  deviations;  vertical 


introduced  in  the  Papua  Region.  That  law  became  uniform  in 
the  New  Guinea  Region  in  1975  (animals  > 2.2  m TL  are 
excluded  in  Tables  3-3  - 3-10) . The  subsequent  increase  in 
the  number  of  breeding  adults  significantly  enhanced 
production  in  that  area.  I assume  the  sudden  sustained 
shift  in  harvest  structure  of  £.  novaeqwineae  reflects  the 
increased  availability  of  young.  A similar  plateau  leap  in 

Province  in  1976.  During  that  year  a new  exporter,  Laki 
Trading,  provided  an  outlet  for  watch  strap'  sited  animals 
that  reduced  the  selective  bias  against  very  small 


To  counter  the  economic  waste  in  trading  of  skins  from 
small  juveniles,  a national  buying  scheme  was  introduced  in 
1979  to  encourage  the  sale  of  small  live  crocodiles  that 
could  be  reared  to  maximum  profitable  size  in  commercial 
facilities  (Whitaker  and  Kemp,  1981) . An  immediate  impact 
of  that  action  was  a minimum  eleven-fold  projected  revenue 
increase  from  live  animals  bought  during  1979  and  1980 
(Hall , unpubl . data) . It  also  reinforced  the  increased 
proportion  of  juveniles  in  the  total  harvest.  Hide 
trafficking  in  animals  less  than  18  cm  BW  (1.0  m TL)  became 
illegal  in  1981  (Parker,  1981) . Falsified  purchase  records 
were  partly  responsible  for  the  huge  increase  of  juvenile  C. 
novaeaulneae  harvest  (and  consequent  depressed  harvest  of 


subadults  and  adults)  in  the  Sepik  during  1981  (F.  Whitam, 

Other  factors  important  in  enhancing  the  economic 
return  of  the  crocodilian  resource  were  extension  activities 
aimed  at  improving  skinning  and  preservation  techniques  and 
the  elimination  of  marketing  less  valuable  hornback  (ventral 
incision)  hides.  The  latter  formerly  occurred  in  high 


Harvest  Mortality  Rates 

There  are  no  implicit  assumptions  regarding  the 
stability  of  population  sizes  in  my  calculations  of  harvest 
mortality  rates  (qa)  for  hunted  crocodile  populations 
(Tables  3-6  and  3-9) , nor  do  I assume  that  the  distribution 
of  the  probabilities  of  age-specific  mortality  are  the  same 
for  harvesting  mortality  as  they  are  for  all  other  sources 
of  mortality  combined.  However,  I believe  that  crocodile 
populations  probably  increased  in  response  to  conservation 
measures  during  the  1969-1980  interval.  This  statement  is 
supported  in  part  by  the  fact  that  the  average  harvested 
skin  size  of  both  species  grew  steadily  larger  during  this 
period  (CHAPTER  2),  despite  an  increasing  proportion  of 

Furthermore,  1 assume  that  commercial  harvesting 
represented  the  primary  (but  not  only)  cause  of  mortality 


for  animals  up  to 
against  harvest  e 


TL  in  hunted  populations  (except 
age  for  which  a selective  bias 


additive  is  unknown  and  unaccounted  for.  Thus,  my  mortality 
estimates  represent  maximum,  not  absolute,  values.  The 
utility  of  such  calculations  lies  in  that  their  allowing 
comparative  relative  estimates  of  the  impact  of  harvesting 
on  populations  within  and  among  species. 

A noted  difference  in  the  harvest  of  these  species  is 
that  the  percentage  of  juvenile  £.  porosus  is  enormously 
higher  than  that  of  £.  nqvaequinsae:  the  juvenile  to 

rest  ratio  is  31.9  to  1 in  the  former  species 


versus  4 . 4 to  1 in  the  latter  species . This  fact  is  due  in 
large  part  to  the  greater  demand  for  hides  of  £.  porosus . 
Skins  of  that  species  produce  a superior  leather  as  a result 
of  their  smaller  belly  scales  and  consequent  greater 
pliability.  Thus,  they  are  worth  about  25%  more  on  the 
domestic  market.  To  a lesser  extent  the  disparate 
juvenile/subadult  harvest  ratio  among  species  is  a 
reflection  of  the  differential  number  of  size  classess 
encompassed  by  such  comparisons.  However,  although  a much 
smaller  proportion  of  £.  porosus  mature  and  breed  below  the 
maximum  legal  harvest  size  compared  with  £.  novaequineae. 
populations  of  £.  porosus  are  exploited  more  intensively 


relative  protection 


low  hunting  mortality  rate  for  £.  porosus  observed  in  New 
Britain  compared  to  mainland  PNG  populations  is  a reflection 
of  the  small  number  of  0-1  year  age  class  animals 
represented  in  the  harvest.  I believe  this  to  be 
attributable  to  marketing  infrastructure  in  that  area, 
although  it  in  part  may  also  reflect  a real  difference  in 
population  size  structure. 

A harvest  feature  common  to  both  these  species  is  the 
disproportionate  number  of  individuals  contained  in  the 
largest  legal  size  class  (see  data  for  20"  BW  row  in  Table 
B-l,  Appendix  B) . Virtually  all  of  these  irregularities 
were  attributable  to  a single  exporter  (Hall,  unpubl.  data) 
and  are  believed  to  be  the  result  of  illegal  trafficking  in 
oversized  hides  carefully  shrunken  to  meet  legal 


PNG  crocodile  populations  of  animals  < 2.2  m TL  were 
perhaps  subjected  to  annual  harvest  rates  of  20  percent  or 
more  during  the  1969-1980  interval  with  no  apparent  long 
term  negative  biological  effects.  Such  rates  may  appear  to 
be  high  and  excessive  to  some  critics  of  crocodilian  trade 
practices.  However,  continuing  research  in  the  southeastern 
U.S.  on  recovered  populations  of  American  alligator 

) compares  favorably  with  the  PNG 


experience.  Studies  designed  to  test  the  effects  of  a 50% 
removal  rate  of  annual  production  (eggs  and  hatchlings)  in 


several  areas  of  northcentral  Florida  indicate  that 
populations  to  date  have  not  been  adversely  affected  by  such 
harvests  (Jennings  et  al.,  1988).  Also,  the  fact  that  PNG 
harvests  were  heavily  directed  toward  juvenile  and  subadult 
animals  accords  well  with  the  life  history  strategies  of 
crocodilians  as  these  are  the  most  expendable  population 

Qfeservefl  Sex  Ratios 

The  sex  of  C.  porosus  and  all  other  crocodilian  species 
studied  to  date  is  known  to  be  temperature-dependent  during 
incubation  (Ferguson  and  Joanen,  1982;  Webb  et  al.,  1963a; 
Webb  and  Smith,  1984;  Lang,  1987).  Presumably  this  is  true 

climate  and  microhabitat  could  be  responsible  for  the 

Middle  Fly  and  Lake  Murray  subdistricts.  Annual 
precipitation  in  the  former  area  (Bosset  = 258  cm;  Hall, 
1983)  is  considerably  lower  than  in  the  latter  (Baboa  = 306 
cm,  Pangoa  - 325  cm;  McAlpine,  1971;  Hall,  1983).  Whether 
possible  differences  also  existed  with  respect  to 
age/sex-specific  mortality  remains  unknown. 

Areas  for_Eurther  Research 

Most  of  the  results  described  herein  are  rather 
straightforward  and  readily  analyzed  using  data  compiled  by 
buyers  during  regulated  business  transactions.  However,  the 


estimation  of  harvest  mortality  rates  (Tables  3-6  and  3-9) 
entails  the  use  of  accessory  data  on  size/age-specific 
growth  rates  and  size/age-specific  sex  ratios  that  are  not 
obtainable  without  special,  and  more  costly,  effort.  In 
order  to  derive  baseline  approximations  of  those  values,  I 
have  had  to  sacrifice  accuracy  for  precision  and  make 
compromising  guesstimates  with  the  best  available  data.  To 
wit,  growth  and  sex  ratio  data  were  either  lacking  or  too 
limited  to  be  useful  for  northern  (New  Guinea  Region)  C. 

Thus,  data  on  growth  and  sex  ratios  for  southern  (Papua 
Region)  £.  novaeau i neae  were  substituted  as  approximations 
for  northern  population  estimates.  Likewise,  growth  rates 
for  PNG  C.  porosus  populations  had  to  be  estimated  by  using 
data  from  Australian  Q.  porosus  populations,  while  sex  ratio 
distributions  for  that  species  were  assumed  to  be  similar  to 

populations  is  also  quite  limited  for  individuals  > 1.5  m 
TL.  Clearly,  additional  investigation  of  intra-  and 
inter-specific  population  growth  rates  and  sex  ratio 
compositions  and  their  associated  variances  is  called  for  to 
improve  the  pictures  presented  herein  and  the  acquisition  of 
such  parameters  should  be  part  of  crocodilian  research 


programs  everywhere. 


foregoing 


The  need  for  greater  accuracy  of  data  for  the 
interpretations  is  best  illustrated  with  respect  to  PNG 
populations  of  £.  novaeauineae.  Although  the  range  of  sizes 
harvested  is  similar  between  northern  and  southern 
populations,  there  are  important  differences  among  them. 
Southern  £.  novaeauineae  are  significantly  larger  at 
hatching  and  hatch  over  a much  longer  interval  than  do 
northern  animals  (Hall  and  Johnson,  1987) . However,  maximum 
verified  sizes  (TL)  of  male  and  female  individuals  are 
smaller  for  southern  £.  novaeauineae  (3.31  m and  2.65  m, 
respectively;  Hall  and  Johnson,  1987)  than  for  their 
northern  counterparts.  A captive  female  at  Pagwi,  East 
Sepik  Province  measured  3.05  m (Jelden,  1985)  while  a record 
sized  male  of  3.47  m was  recently  found  illegally  killed  in 
the  Mamberamo  River,  Irian  Jaya  (the  western,  Indonesian, 
half  of  New  Guinea)  (S.  Frazier,  cited  in  Brazaitis  et  al., 
1988) . Additionally,  and  more  importantly,  subsequent  to 
the  compilation  of  the  foregoing  harvest  data,  it  was  found 
that  the  southern  population  of  £.  novaeauineae  was 
morphologically  distinct  from  the  northern  type  population 
(Hall,  1989) . 
summary 

Although  the  PNG  crocodile  harvest  has  been  sustained 
with  no  detectable  long-term  negative  effects,  the  potential 
for  further  economic  return  from  a manipulated  harvest  is 


Johnson  ( 1987 ) 


female  c . novaeouineae  In  Western  Province  were  smaller  than 
the  size  protected  by  law.  Accordingly,  they  recommended  a 
reduction  in  the  upper  skin  size  limit  to  40  cm  BW  (<  1.8  m 
TL)  to  enhance  juvenile  production.  Cox  (1985}  also 
reported  that  significant  proportions  of  breeding  female  c. 
novaeouineae  (86%)  and  that  some  £.  porosus  were  unprotected 
in  the  Sepik  Provinces  as  well.  Consequently,  he  made  a 
similar  suggestion  to  reduce  the  upper  size  trading  limit  as 
did  Montague  (1983). 

In  conclusion,  differential  harvesting  and  marketing 
strategies  have  probably  acted  upon  regional  populations  of 
PNG  crocodile  species  to  make  them  appear  more  dissimilar 
than  what  they  really  are.  These  factors  can  largely 
account  for  trends  detected  in  size  structure  within  £. 
novaeouineae  harvests  and  for  differing  intra-  and 
inter-specific  harvest  mortality  rates.  Because  of  the  very 
large  sample  sizes,  the  significant  differences  observed 
among  intra-specific  comparisons  of  harvest  compositions  are 
less  meaningful  biologically  than  they  are  statistically. 
This  underscores  the  need  for  caution  when  interpreting  this 

The  people  of  PNG  have  had  a long  cultural  affinity 
with  crocodiles.  That  fact,  coupled  with  large  tracts  of 
difficult  access  habitat,  low  density  human  populations,  and 


CHAPTER  4 

ESTIMATION  OF  NESTING  FEMALE  CROCODILIAN  SIZE  FROM  CLUTCH 
CHARACTERISTICS.  CORRELATES  OF  REPRODUCTIVE  MODE  AND 
HARVEST  IMPLICATIONS 


(Sill, 


latirostris . C.  vacare)  (Neill,  1971;  Groombridge,  1982) . 


The  surviving  members  of  the  genus  Alligator  (a* 
Blaaigg.lPPiM8.l8  and  a.  sinensis)  occur  only  in  temperate 


Guinea  crocodiles  (C.  novaeguineae) . The  former  species  is 


Guinea  central  highlands,  as  they  compose  a taxonoraically 
distinctive  form  (Hall,  1989) . Biologically  meaningful  data 
on  reproductive  female  crocodilians  are  generally  difficult 
to  obtain  compared  to  data  on  crocodilian  nests  and  nest 
sites.  Therefore,  I sought  to  determine  whether  the  latter 
type  of  data  could  be  used  to  reconstruct  demographic 
profiles  of  breeding  crocodilian  populations  given  data  on 
small  subsets  of  the  nesting  female  segments.  I also 
briefly  consider  the  implications  of  harvesting  on  the 
reproductive  demography  of  these  populations. 

study  Areas 

in  Lake  Murray  District,  Western  Province,  Papua  New  Guinea 
(PNG;  7‘00'S,  141'30'E)  (Hall  and  Johnson,  1987;  Hall, 
unpubl.).  Additional  data  on  nesting  animals  are  from  Cox 
(1985:  Appendix  6).  The  Lake  Murray  area  is  primarily 
tropical  moist  forest  (mean  annual  precipitation  3,200  mm) 
with  pronounced  wet  (November-May)  and  dry  (June-October) 
seasons.  Habitats  important  for  nesting  crocodiles  are  the 
fringing  mixed  swamp/woodland  forests,  Echinochloa  grass 
swamps,  and  Saccharum-Phragmites  grass  swamps.  Further 
details  are  provided  in  CSIRO  (1971)  and  Hall  and  Johnson 

Observations  on  nesting  alligators  were  made  at 
Lacassine  National  Wildlife  Refuge  (NWR;  30'00'N,  92'50'W), 


Louisiana  in  1985  and  . 

»t  Lake  Griffin  (28‘50'N,  81‘50'W), 

Lake  Okeechobee  (27 ‘00 

'N,  81  * 00 1 W) , Lake  Woodruff  NWR 

(29-05'N,  81'25'W) , an. 

i orange  Lake  (29’30'N,  82‘10'W), 

Florida  during  1987.  Lacassine  NWR  (12,860  ha)  is  within 
the  freshwater  coastal  marsh  region  of  Louisiana  and  about 
one-half  of  the  study  site  was  a shallow  (1-2  m)  impoundment 
containing  large  stands  of  maidencane  (Panicum  hemitomonl 
and  bulltongue  tsaaittaria  falcatal  that  supported  dense 


nesting  by  alligators 

(an  estimated  1 nest  per  13  ha) . 

Fruge  (1974)  provides  i 

i full  synopsis  of  the  area's 

vegetation.  Descriptions  of  the  Florida  study  areas  are 
contained  in  Deitz  and  Hines  (1980;  Orange  Lake),  Jennings 
et  al.  (1987;  Lake  Griffin,  Lake  Woodruff  NWR) , and 
Millerson  (1987;  Lake  Okeechobee). 


Methods 

The  morphometric  < 

»nd  reproductive  variables  recorded, 

their  nearest  unit  of  i 

neasure,  and  their  abbreviations,  are 

as  follows:  body  mass  (0.1  kg)  = BM,  clutch  mass  (g)  = CM, 
clutch  size  ■ CS,  dorsal  cranial  length  (0.5  cm)  [measured 
from  the  anterior  snout  to  posterior  supraoccipital ] ■ DCL, 
snout-vent  length  (cm)  = SVL,  maximum  tail  girth  (cm)  ■*  TG, 


total  length  (cm)  - T03 

'.  SVL  measurements  for  crocodiles 

were  taken  from  the  snc 

>ut  tip  to  the  anterior  margin  of  the 

cloaca  as  per  Australia 

in  convention  (Webb  and  Messel,  1978), 

to  the  posterior 


alligators 


recorded 


cloaca  consistent  with  the  North  American  standard  (Chabreck 
and  Joanen,  1979) . Morphometric  measurements  of  alligators 
were  from  animals  captured  by  short  penetration  harpoon 
(modified  design  of  Webb  and  Messel,  1977a)  or  from 
individuals  live-trapped  at  nest  sites  via  an  unbaited 
peddle  trip  snare  (Mazzotti  and  Brandt,  1988) . 

Additionally,  I examined  reproductive  tracts  of  female 
alligators  harvested  at  Lacassine  NWR  in  August-September 
1985  to  compare  sizes  of  active  breeding  animals  with  those 
estimated  from  nest  studies.  Females  were  classified  as 
sexually  mature  based  on  the  color  and  texture  of  their 
ovaries  (Graham,  1968;  Joanen  and  McNease,  1980). 
Reproductive  activity  was  assessed  by  the  presence  or 
absence  of  postovulatory  ovarian  follicles  (Lance,  1989) . 

Nesting  crocodiles  were  not  captured,  although  their 
clutch  parameters  were  recorded.  Morphometric  measurements 
were  estimated  as  follows.  Data  from  Cox  (1985,  Appendix  6) 
were  supplemented  to  that  of  Hall  and  Johnson  (1987)  to 
derive  a least  squares  relationship  between  CM  and  C5  for 

estimate  CMs  where  unknown  for  individual  clutches  (N  = 24). 
X also  regressed  CS  against  nesting  female  TOT  (modified 
from  Hall  and  Johnson,  1987)  to  estimate  female  TOTs  for 
where  the  adult  was  unobserved  (N  = 55),  whereby: 


clutches 


Estimated  TOTs  were  then  transformed  to  SVL  estimates  by 
using  equation  2 of  Montague  (1984b):  SVL  ■ 5.73  + 0.49  TOT 

The  method  for  analysis  of  £.  novaeouineae  oviposition 
season  data  is  described  by  Caughley  (1977:72-78).  The 
dates  of  nest  detection  were  corrected  to  the  dates  of  egg 
deposition  through  the  use  of  embryo  growth  curves  (Hall. 


Alligator  m 


ssippisnsls 


The  morphometric  data  on  female  alligators  from  Florida 
associated  with  known  clutch  data  (N  - 11)  were  too  few  to 
permit  meaningful  comparisons  among  study  sites;  however, 
ween  areas  were  similar  to  each  other 
s pooled  for  analyses.  Since  among 
e strong,  equations  derived  from  those 
o calculate  snout-vent  length  (SVL) 

3 Louisiana  animals  for  which  only 

was  chosen  as  a more  relevant  indicator  of  adult  female  size 
than  total  length  (TOT)  due  to  the  high  proportion  of  tail 
injuries  among  crocodilians  (Webb  and  Messel,  1977b; 
Montague,  1984a;  Magnusson,  1985). 

The  best  variable  for  predicting  the  size  of  an  adult 
from  a single  attribute  e 


regression  slopes  b 


Florida  a 


estimates  ol 


female  alligator 


animal  was  CM,  although  the  dorsal  cranial  length  (DCL) 
measurement  of  an  isolated  skull  was  a slightly  more 
sensitive  predictor  of  SVL  (Table  4-1,  Fig.  4-1) . The 
relations  between  these  and  other  tested  factors  were 
essentially  isometric  and  were  not  significantly  enhanced  by 
the  use  of  logarithmic  transformations. 

Allometric  growth  in  skull  length  to  body  proportions 
within  the  Crocodylia  has  been  documented  (Kalin,  1933;  Webb 
and  Messel,  1978;  Hall,  198Sa) , but  this  was  not  observed 
within  the  size  range  (97-140  cm  SVL)  of  female  alligators 
in  this  study.  A trend  toward  allometry  of  DCL  would 
probably  be  present  in  larger  animals,  but  may  be 
sex-specific  and  outside  the  body  lengths  attained  by  female 
alligators. 

Where  morphometric  data  on  the  associated  females  were 
lacking,  I derived  SVLs  for  nesting  female  alligators  from 
Lake  Woodruff  NWR,  orange  Lake,  and  Lacassine  NWR  from  their 
CMs  (equation  3,  Table  1).  Estimated  sizes  (x  = 119.6  cm) 
did  not  differ  (t  - -6.37,  R - 11,  P > 0.05)  for  animals  on 
which  actual  sizes  (x  = 119.5  cm)  were  obtained.  Thus, 
estimates  were  pooled  with  captured  nesting  female  data  to 
develop  demographic  profiles  of  breeding  females  from  these 
populations  (Fig.  4-2) . 

Nesting  alligators  (pooled  values  of  actual  and 
estimated  sizes)  at  Orange  Lake  (x  SVL  = 106.0  + 7.39  SD  N 


slightly  smaller 


CLUTCH  MASS  (g) 


Figure  4-1.  Regression  plot  of  clutch  mass  and  nesting 
female  alligator  snout-vent  length. 


SNOUT-VENT  LENGTH  (cm) 


Figure  4-2.  Size  structures  of  nesting  female  Alligator 
mississippiensis  from  harvested  (Orange  Lake  and  Lacassine 
National  wildlife  Refuge  - top)  and  protected  (Lake  Woodruff 
National  Wildlife  Refuge  - bottom)  populations. 


> 0.05),  animals  at  Lacassine  NWR  (x  = 110.1  + 8.89,  N = 9) . 

harvest.  However,  reproductive  females  at  Orange  Lake  (t  - 
3.01)  and  Lacassine  NWR  (t  ■ 2.93)  were  significantly 
smaller  (P  < 0.01)  than  alligators  nesting  at  Lake  Woodruff 
NWR  (x  = 120.1  ± 7.49,  N = 18) , which  was  a protected 
population.  Furthermore,  the  sizes  of  nesting  alligators 
from  Lacassine  NWR  did  not  differ  (t  ■ 1.43,  P > 0.05)  from 
that  of  harvested  females  (x  - 104.0  + 11.95  SD,  N ■ 40)  who 
had  nested  that  year,  based  on  reproductive  tract  analysis. 
Thus,  those  size  distributions  were  included  in  Figure  4-2. 
Although  sample  sizes  from  Orange  Lake  are  small,  those 
estimates  are  consistent  with  observations  on  the  mean  size 
of  sexually  mature  female  alligators  removed  from  those 
systems  during  regulated  harvest  (A.  Woodward,  pers.  comm. ) . 
The  Lake  Woodruff  NWR  sample  encompassed  about  one  third  of 
the  known  nesting  segment  for  that  season  so  those  estimates 
are  probably  not  strongly  biased. 

X converted  estimated  SVL  measurements  for  Louisiana 
alligators  to  TOT  estimates  by  using  equation  14  of  Table 
4-1  which  was  developed  from  Florida  data.  The  same  SVL 
values  were  then  compared  with  transformed  TOT  estimates 
derived  from  Louisiana  observations,  whereby  female  SVL  = 
-2.24  + 0.5035  TOT  and  body  mass  (BM)  = 1.35  - 0.0378  TOT  + 
0.0000046  TOT3  (chabreck  and  Joanen,  1979).  Those  results 
indicated  that  female  alligators  from  Louisiana  possess 


longer  TOTs  than  do  Florida  animals  of  an  equivalent  SVL; 
hence,  with  corresponding  SVL,  alligators  from  Louisiana  had 
a heavier  BM  than  Florida  specimens,  consequently,  the 
relative  clutch  mass  (expressed  as  the  ratio  of  ([(CM  / 

1000)  / BM]  x 100)  was  depressed  in  Louisiana  females  (x  - 
5.12  + 0.41  SD,  range  - 4.7-5.B)  when  contrasted  with 
Florida  alligators  (x  ■ 6.82  + 1.14,  range  a 5. 1-8. 7).  This 
difference  (t  = 4.01,  P < 0.001)  might  seem  to  be  an 
artifact  of  the  larger  mean  SVL  (119  cm)  of  females 
composing  the  Florida  sample.  However,  the  comparable 
relative  clutch  mass  for  a Louisiana  individual  with  an  SVL 
of  121  cm  was  only  5.6,  suggesting  that  the  difference  is 


The  estimated  mean  SVL  for  60  individuals  (excluding 
four  captives)  was  104.1  cm  ± 10.89  SD  (Fig.  4-3). 

Estimated  CM  and  BM  (equation  83  of  Montague,  1984b)  for  a 
crocodile  of  this  size  were  2,001  g and  38.2  kg, 
respectively.  Relative  clutch  mass  was  5.2,  which  lies 
within  the  interval  of  4-7%  previously  reported  for  the 
species  (Hall  and  Johnson,  1987). 

The  lower  mean  SVL  of  C . novaequlneae  compared  to 
alligators  was  a reflection  of  smaller  adult  size.  Maximum 
recorded  size  (TOT)  for  female  £.  novaeonineae  is  2.65  m 


Figure  4-3.  Size  structure  of  nesting  female  Crocodvl us 
novaequineae  in  Lake  Hurray  District,  Papua  New  Guinea. 


(Hall  and  Johnson,  1987),  whereas  female  alligators  can 
attain  2.87  m (Joanen  and  McNease,  1989). 

The  sizes  (SVL)  of  nesting  crocodiles  reported  here 
encompassed  the  minimum  and  maximum  known  range  of 
reproduction  by  the  species.  Since  data  representing  the 
upper,  but  not  the  lower,  limit  of  the  reproductive  curve 
were  available  for  alligators,  I hypothesized  that  an 
estimate  of  the  minimum  size  at  reproduction  in  alligators 
would  be  proportional  to  the  maximum  adult  female  size  ratio 
between  these  species  (0.92) . To  test  this  assumption,  I 
took  the  lowest  observed  SVL  value  for  £.  novaequlneae  (80 
cm) , added  3 cm  to  make  the  measurement  equivalent  to  that 
for  alligators  (see  METHODS) , and  obtained  a predicted  SVL 
estimate  of  90  cm  (83  / 0.92)  for  alligators.  This  value 
closely  approximated  that  observed  (88  cm)  for  the  smallest 
reproductively  active  female  alligator  that  I found  during 
reproductive  tract  analysis  of  harvested  animals  from 
Lacassine  NWR  and  is  similar  to  that  for  a 1.71  m TOT  female 
observed  during  incubation  (Joanen  and  McNease,  1989) . 

Thus,  in  these  two  crocodilians,  females  appear  capable  of 
nesting  once  they  have  attained  about  55-60*  of  maximum  TOT. 

Existing  legislation  designed  to  safeguard  breeding 
population  segments  protects  PNG  crocodiles  above  51  cm 
commercial  belly  width  from  harvest.  This  measurement  is 
the  equivalent  of  an  animal  with  an  SVL  of  108  cm  (209  cm 
TOT)  (equation  12  of  Montague,  1984b) . However,  based  on 


the  analysis  presented  here,  57%  (95%  confidence  intervals  - 
40-73%)  of  nesting  £.  novaeguineae  were  below  that  size.  A 
previous  estimate  of  74%  (Hall  and  Johnson,  1987)  was  based 
on  visual  estimates  of  attendant  nesting  females.  The 
discrepancy  among  these  estimates  is  due  to  the  clustering 
around  the  minimum  size  at  protection  (see  Fig.  4-3) . 
Elsewhere,  Cox  (1985)  found  that  86%  of  female  £. 

w Guinea  killed  at  their 
e protected  by  law. 


The  nesting  season  for  £.  novaeguineae  from  Lake  Murray 
is  almost  exclusively  coincidental  with  the  wet  season 
(Hall,  1985b) . Unlike  American  alligator  populations  that 
confine  egg  laying  activities  to  a period  of  approximately 
three  weeks  during  June  or  early  July  (Joanen  and  McNease, 
1979) , Lake  Murray  crocodiles  oviposit  during  an  interval 
more  than  ten  times  as  long  (Fig.  4-4) . The  mean 
oviposition  date  was  February  6 ± 345  days  SD,  far  in  excess 
of  a 30  day  SD  characteristic  of  populations  that  follow 
birth  flow  models  (Dittus,  1981) . Clutch  size  remained 
nearly  constant  throughout  the  oviposition  period. 


This  study  clearly  demonstrates  that  demographic 
interpretations  of  the  nesting  female  segments  of 
crocodilian  populations  are  possible  on  the  basis  of 
associated  clutch  parameters.  Although  nesting  densities 


crocodilian  populations 


(McNease  and  Joanen,  1978;  Hollands,  1987),  monitoring  of 
the  reproductive  cohorts  per  se  has  not  been  possible  until 
now  due  to  the  paucity  of  data  on  breeding  females 
associated  with  specific  nests.  The  methods  I used 
circumvented  that  difficulty  and  permitted  a reliably 
accurate  determination  of  the  size  of  nesting  females  among 
populations  from  existing  data  sets  composed  of  CM  or  CS 
alone.  These  results  indicate  that  CM  is  a better  predictor 
of  adult  female  size  than  CS  for  both  alligators  and 
crocodiles.  They  lend  further  support  to  the  use  of  SVL 
measurement  as  the  most  biologically  meaningful  criterion  of 
crocodilian  size  since  female  alligators  with  equivalent 
SVLs  from  geographically  disjunct  populations  are  not 
equivalent  in  TOTs.  This  phenomenon  was  also  reported  for 
Crocodvlus  porosus  in  which  Bornean  crocodiles  possessed 
longer  tails  than  did  Australian  animals  (Webb  and  Messel, 

Temperate  alligators  and  tropical  crocodiles  share 
similar  behavioral  repertoires  (Garrick  and  Lang,  1977) . 
However,  the  timing  and  duration  of  reproductive  seasonality 
differs  strongly  among  fi.  mississjppiensis  and  southern-form 
£.  novaequineae. 

The  full  complement  of  reproductive  activities 
(courtship  through  hatching)  occurs  between  mid-March  and 
mid-September  in  fi.  mississjppiensis  (Ferguson,  1985), 


whereas  reproduction  in  £.  novaequineae  is  circa-annual 
(Hall  and  Johnson,  1987).  In  addition,  subsets  of  the 
reproductive  cycle  are  temporally  compressed  in  alligators — 
incubation  requires  9-10  weeks  to  completion  (Joanen,  1969) 
versus  12-13  weeks  for  New  Guinea  crocodiles  (Hall  and 
Johnson,  1987) . Consequently,  the  genitalia  of  hatchling 
alligators  are  morphologically  undifferentiated  compared  to 
hatchling  crocodiles.  Thus,  neonate  alligators  cannot  be 
reliably  sexed  by  cloacal  examination  (Joanen  and  McNease, 
1978) , although  gender  determination  of  hatchling  crocodiles 
by  that  method  is  possible  (Webb  et  al.,  1984).  clearly, 
the  narrow  time  interval  for  egg  laying  within  alligator 
populations  characterizes  the  species  as  having  a pulse 
reproductive  mode  (caughley,  1977) , contrasted  with  the  flow 
model  of  £.  novaequineae.  Why  these  species  differ  in  this 
respect,  and  what  some  of  the  probable  consequences  of  each 
mode  on  population  recruitment  are,  remain  to  be  answered. 

The  timing  of  oviposition  by  American  alligators  is 
negatively  correlated  with  ambient  temperatures  during  the 
spring  months  of  March  to  May  (P  < 0.01;  i.e.,  years  of 
colder  mean  spring  temperatures  result  in  later  mean  dates 
of  egg  laying) . Egg  deposition  also  coincides  with  the 
maximal  photoperiod  interval.  Although  rainfall  does  not 
affect  the  time  of  nesting  it  does  influence  the  degree  of 
nesting  (Joanen  and  McNease,  1979;  Kushlan  and  Jaoobsen, 

1990) . This  agrees  with  what  one  should  expect  of  a 


poikilothermic  species  in  temperate  habitat  where  marked 
seasonal  fluctuations  in  incident  sunlight  and  temperature 
greatly  constrict  the  period  of  primary  and  secondary 
productivity  compared  to  tropical  environments.  It  is  also 
not  surprising  that  alligators  have  a size-specific 
dominance  breeding  hierarchy  whereby  larger  males  and 
females  copulate  earlier  and  more  frequently  within  a season 
than  do  smaller  conspecifics  (Ferguson,  1965;  Vliet,  1987) , 
given  the  limited  temporal  framework  within  which  mating 
must  occur  (two  months  annually) . Accordingly,  model 
simulations  for  alligator  populations  are  largely 
deterministic  since  the  proximate  factors  affecting 
recruitment  are  essentially  stable  and  predictable,  whereas 
stochastic  events  (such  as  hurricanes)  are  relatively  rare 
(Nichols  et  al.  1976).  However,  the  impact  of  stochasticity 
on  alligator  populations  can  be  devastating.  For  instance, 
the  early  arrival  of  1985  hurricanes  in  Louisiana  resulted 
in  the  loss  of  about  80-90%  of  the  coastal  alligator  egg 
crop  due  to  flooding  (T.  Joanen,  pers.  comm.).  The 
long-term  effects  of  such  occurrences  on  the  population  are 
undoubtedly  lessened  because  the  potential  longevity  of 
alligators  (and  other  crocodilians)  may  exceed  50  years 
(Chabreck  and  Joanen,  1979,  Graham,  1968).  Thus,  breeding 
segments  of  crocodilian  populations  may  be  composed  of  20  or 
more  active  age  cohorts.  Some  female  crocodilians,  though, 
eventually  become  reproductively  senescent 


(Graham, 


Joanen  and  McNease,  1980) . For  further  discussion  on  the 
Abercrombie  (1989) . 

In  contrast,  the  proximate  factor  apparently  governing 
the  reproductive  cycle  of  £.  novaeouineae  is  annual 
precipitation,  as  photoperiod  and  temperature  are  fairly 
uniform  throughout  the  year.  Pronounced  variation  in 
rainfall  occurs  both  within  and  among  years,  with  most  of 
the  variance  accounted  for  by  patterns  of  dry  season 
precipitation  (Hall,  chapter  3).  Although  the  reproductive 
hormonal  cycle  has  yet  to  be  studied  in  this  species, 
courtship  behavior  appears  to  be  synchronous  with  the  onset 
of  the  wet  season  rains  as  egg  deposition  spans  the  duration 
of  the  wet  season.  Crocodiles  seem  to  cue  on  the  high  water 
mark  of  the  previous  year  as  their  nests  are  constructed  on 
land  above  that  point  or  are  built  on  mounds  of  floating 
herbaceous  vegetation  that  rise  and  fall  with  water  levels 
(Hall  and  Johnson,  1987) . The  prolonged  seasonality  of 
nesting  makes  it  unlikely  that  catastrophic  stochastic 
events  (e.g.,  earthquakes,  fires,  storms)  will  eliminate  an 
entire  cohort  as  can  happen  in  alligator  populations 
(above) . Whether  or  not  a size-specific  dominance  breeding 
hierarchy  exists  in  wild  £.  novaeguineae  is  also  unknown, 
although  such  behavior  is  exhibited  by  captive  animals 
(Lang,  pers.  comm.;  Hall,  pers.  obs.).  I suspect  that  if 
that  behavior  occurs  in  the  wild,  the  intensity  of 


interactions  should  be  relaxed  compared  to  alligators 


the  interval  in  which  spermatogenesis  and  vitellogenesis  c 
occur  is  potentially  four  times  as  long  i 

Although  techniques  for  aging  crocodilians  an 
long-lived  reptiles  have  been  developed  (Graham,  1 
Hutton,  1986;  Zug  et  al.,  1986;  Zug  and  Rand,  1987 
produce  estimates  that  suffer  from  low  precision  a 
accuracy  f 


t very  young  cohorts.  It  is  partly  for 
this  reason  that  Nichols  (1987)  argued  for  the  use  of  size 
rather  than  age-specific,  cohort  analysis  of  crocodilian 
populations.  In  addition,  a population  with  a flow  mode  o 
reproduction  like  southern-form  £.  novaeauineae  would 
further  obfuscate  such  interpretations,  especially  if  maki 
inter-population  comparisons  with  northern-form  £. 

las  a much  shortened  reproductive  seaso 


(Cox,  1985) . Since  reproductive  mode  may  differ  among 
populations  within  a species,  these  generalizations  should 
not  be  assumed  to  uniformly  apply  to  other  crocodilians  on 
the  basis  of  latitudinal  gradients. 

I believe  that  the  effects  of  controlled  harvests  may 
have  operated  differentially  on  the  breeding  female  segments 
of  these  populations  of  4.  mississippiensls  and  £. 
npyaequlneae-  The  similarity  of  mean  sizes  of  nesting 
female  alligators  in  populations  exposed  to  harvest 
(Lacassine  NWR  and  Orange  Lake)  and  the  significantly  larger 
mean  size  of  animals  nesting  in  an  unhunted  population  (Lake 


Woodruff  (TOR)  suggests:  (X)  the  existence  of  a reservoir  of 
sexually  mature,  but  socially  subordinate,  animals  in 
recovered,  protected  alligator  populations;  (2)  that  such 
animals  fail  to  nest  because  they  are  exploitatively  or 
behaviorally  inferior  to  larger,  dominant  females;  and,  (3) 
that  they  will  not  achieve  active  reproductive  status  until 
either  attaining  larger  sire  or  until  a socially  mature 
female  is  removed  through  natural  mortality  or  managed 
harvest,  thus  producing  a reproductive  vacancy.  The  mean 
size  of  nesting  females  in  the  protected  population  was  14% 
longer,  but  72%  heavier,  than  that  of  females  in  harvested 
populations.  Thus,  even  relatively  small  differences  in 
length  can  place  smaller  sexually  mature  females  at  a great 
competitive  disadvantage  with  larger  conspecifics. 

A caveat  of  the  preceding  paragraph  is  the  implicit 
assumption  that  competition  among  females  is  occurring  due 
to  some  strongly  limiting  factor,  presumably  access  to 
nesting  habitat.  Habitat  quality  differs  greatly  throughout 
the  range  of  these  species,  where  nesting  densities  are 
generally  low  in  riverine  areas,  intermediate  in  lake 
habitats,  and  highest  in  grassland  marshes.  Among  the 
latter , Eanicum  hemitomom  marsh  is  considered  qualitatively 
inferior  to  coastal  brackish  and  intermediate  marshes  in 
Louisiana  (T.  Joanen,  pers.  comm.).  Nesting  habitat  appears 
not  to  be  a limiting  factor  in  these  marshes  and  they 
support  the  highest  nesting  density  of  alligators  in  the 


United  States.  The  percent  of  annual  nesting  by  adult 
female  alligators  in  these  marshes  is  over  twice  as  high 
(68%;  Joanen  and  McNease,  1980J  as  reported  in  Florida  (29%, 
Kushlan  and  Jacobsen,  1990)  or  South  Carolina  (27.5%, 
Wilkinson,  198%).  Accordingly,  management  harvest  quotas  in 
Louisiana  are  based  on  habitat  type  (Joanen  and  McNease, 

These  results  also  suggest  that  selective  culling  of 
larger  nesting  female  alligators  may  lead  to  accelerated 
population  growth  through  the  creation  of  higher  nesting 
densities  and  incremented  egg  and  hatchling  production  (see 
Hines  and  Abercrombie,  1987  and  Aberccrombie,  1989) . clutch 

(Ferguson,  1985)  and  £.  novaeguineae  (Hall  and  Johnson, 

1987) . A decrease  in  mean  clutch  size  resulting  from  the 
prudent  removal  of  large  nesting  female  alligators  should  be 

absence  of  more  dominant  females. 

The  relatively  short,  clearly  demarcated  nesting  season 
of  alligators  (in  addition  to  differential  seasonal  habitat 
use  by  adult  males  and  females)  favors  management  practices 
that  selectively  cull  adult  animals  following  the  completion 
of  reproduction.  Such  measures  would  not  be  feasible  for 
the  management  of  southern-form  C . novaeguineae  due  to  its 
very  extended  reproductive  period.  Harvest  regulations 
target  alligators  larger  than  1.2  m TOT,  whereas  harvests  of 


crocodiles  are  restricted  to  animals  less  than  2.1  m. 
Although  the  life  history  traits  of  crocodilian  species 
generally  support  the  removal  of  juvenile  and  subadult  size 
classes  (Messel  et.  al,  1981!  Magnusson,  1986),  most 
reproductively  active  £.  novaeguineae  were  nesting  at  sizes 
smaller  than  that  protected  by  law  and  were  thus  also 
exposed  to  harvest.  Consequently,  the  estimated  annual 
turnover  in  the  female  breeding  segment  was  high  (>  60%)  and 
population  growth  retarded  (Hall  and  Johnson,  1987) . 

Additional  studies  of  this  nature  are  needed  to  more 
fully  evaluate  whether  or  not  site-specific  differences 


herein  and  to  assess  their  predictive  ability  for 
constructing  demographic  profiles  of  breeding  female 
crocodilian  populations.  Although  based  on  small  sample 
sizes,  these  models  are  practical  because  they  provide  a 
fair  to  excellent  accounting  of  the  associated  variances  (r2 
values)  for  each  model.  Preliminary  tests  indicate  that 
they  provide  excellent  agreement  of  estimated  with  actual 
sizes  of  nesting  female  alligators  subsequently  captured  (H. 
P.  Percival,  unpubl.  data).  The  use  of  multiple  regression 
analysis  did  not  enhance  those  values  since  predictors  were 
highly  correlated  with  each  other. 


In  conclusion,  accurate  and  precise  estimates  of 
nesting  adult  female  crocodilians  are  feasible  in  the  course 
of  activities  normally  conducted  by  ecologists  and 


management  authorities,  with  little  expenditure  and 
justifiable  effort.  This  has  proven  effective  with  both 
hole  nesting  Crocodvlus  acutus  (Hazzotti  and  Brandt,  1988) 
and  mound  nesting  species  (this  study) . Data  of  these  type, 
hitherto  rarely  obtained,  are  of  demonstrable  value  to 
biologists  and  agency  personnel  concerned  with  crocodilian 
conservation.  An  initial  understanding  of  the  sizes  of 
reproductive  females  provides  a mechanism  for  the  long-term 
monitoring  of  crocodilian  populations  in  response  to 
protection  or  induced  harvest  strategies,  as  well  as 
affording  a practical  manner  in  which  unanswered  questions 
pertaining  to  the  reproductive  longevity  and  success  of 
individual  animals  may  be  addressed. 


CHAPTER  5 

CRANIAL  MORPHOMETRY  OF  NEW  GUINEA  CROCODILES 
: ONTOGENETIC  V 
N RELATIVE  GROWTH  OF  THE  SKULL 
SESSMENT  OF  ITS  UTILITY  A 
THE  SEX  AND  SIZE  OF  INDIVIDUALS 


Introduction 

crocodile  fcrocodvlus  novaeauineae) , endemic  to  the  main 


crocodilians  or  their  byproducts  is  regulated  by  the 


Fauna  and  Flora  (CITES)  (King,  1973)  since  all  23  species  of 
the  order  Crocodylia  (alligators,  caimans,  crocodiles. 
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1980;  Whitaker  and  Kemp,  1981;  Hall,  CHAPTERS  2 and  3) . Due 
to  the  comparative  abundance  of  this  species  and  the  limited 
availability  of  C.  oorosus  material,  I chose  £.  novaecuineae 


Harvest  data  of  PNG  crocodiles  provide  an  indirect 
measure  of  size  (commercial  belly  width;  see  METHODS) , but 
no  information  on  the  sex  or  age  of  individuals.  Gathering 
data  from  freshly  killed  specimens  is  not  feasible  due  to 
rapid  deterioration  of  hides  and  logistic  constraints 
(remote  areas,  manpower,  and  economics).  These  facts 
e funding  of  long-term  studies  to  generate 
data  on  caught  wild  individuals  to  monitor  growth 
d population  sex  ratios.  Such  costs  have  been 


prohibitive  to  all  but  a few  external  research  agencies.  In 
addition,  crocodilians  are  extremely  wary  animals  that 
develop  rapid  learning  avoidance  to  capture  (Bustard,  1968) . 

Consequently,  I sought  to  develop  a cost-effective, 
alternative  method  to  address  these  problems  through 
morphometric  sampling  of  discarded  crocodilian  skulls  from 
hunter  harvests  as  a practical  management  tool  to  aid  in 
monitoring  resource  use.  The  objectives  of  this  study  were 


to:  (1)  examine  and  quantify  ontogenetic  patterns  in 
relative  growth  of  the  skull  in  a crocodilian,  Crocodvlus 
novaeouineae:  and  (2)  assess  the  utility  of  the  functional 
patterns  of  relative  growth  as  a predictor  of  gender.  Such 


knowledge  will  permit  determination  of  population  structural 


parameters  of  future  harvests.  Thus,  demographic  trends  may 
be  monitored  and  the  potential  to  manage  the  resource 
enhanced. 


Background 

Growth  in  the  Order  Crocodylia  is  thought  to  be 
indeterminate  (Cott,  1961;  Jacobsen  and  Kushlan,  1989; 
Brisbin,  1990} . Crocodilians  exhibit  typical  size 
dimorphism,  whereby  males  grow  longer  and  heavier  than 
females  (Neill,  1971) . However,  this  discrepancy  may  not  be 
apparent  where  regulated  crocodilian  harvests  occur  if  the 
upper  size  limit  is  below  the  point  where  the  divergence 
between  sexes  is  pronounced  (such  as  in  New  Guinea) . Other 
than  size,  sexual  dimorphism  in  crocodilians  is  evident  only 
in  socially  mature,  adult  male  gharials  fGavlal is 
oanoetlcus)  that  exhibit  a secondary  bulbous  snout 
protuberance  (Basu,  1980) . However,  the  reliability  of  this 
character  remains  uncertain  as  a necropsied  adult  possessing 
such  a feature  was  found  to  be  female  (Martin  and  Bellairs, 
1977) . Thus,  to  date,  manual  probing  of  the  cloaca 
(Chabreck,  1963)  or  vivisection  (Magnusson  and  Hero,  1990) 
are  the  only  certain  means  of  determining  gender  in  a living 
individual,  except  where  size  dimorphism  is  apparent. 

As  a group,  crocodilians  are  long-lived  (40+  years) 
vertebrates  and  dimorphic  attributes  of  growth  appear  to 


chabreck  and  Joanen,  1979) . The  disparity  of  such  growth  is 
strongly  influenced  by  incubation  temperatures  during 
embryogenesis  (Ferguson  and  Joanen,  1983;  Webb  and  Smith, 
1984) . Functional  changes  in  the  rates  and  directions  of 
growth  for  areas  of  the  crocodilian  skull  are  also  known  to 
occur  with  the  age  and  site  of  individuals  (Dodson,  1975) . 

In  light  of  the  foregoing,  I hypothesized  that  the 
skulls  of  crocodilians  (beyond  the  early  age  classes)  should 
also  exhibit  a dimorphic  suite  of  characters  to  accompany 
the  ontogenetic  changes  in  growth  of  body  size  with  respect 
to  sex.  In  essence,  the  skull  of  a female  crocodile 
compared  to  that  of  an  equal  length  male  of  age  X will 
belong  to  an  animal  of  age  X+t.  Thus,  the  shifts  in 
ontogenetic  patterns  of  relative  growth  of  the  skull  of 
comparable  size  individuals  should  be  more  pronounced  in 
females  than  males  of  equivalent  size,  and  become  even  more 
disparate  as  size/age  increases.  As  such,  the  gender  of 
like-size  individual  skulls  is  indiscernible  with  specimens 
in  hand  (Fig.  5-2) , but  may  become  apparent  when  patterns  of 
growth  are  viewed  and  analyzed  in  a multidimensional 

The  quantitative  study  of  growth  in  crocodilians  is 
fairly  new.  Early  descriptive  investigations  on  aspects  on 
embryonic  growth  include  those  by  Parker  (1883) , Clarke 
(1891),  Voeltzkow  (1902),  Reese  (1908),  Mcllhenny  (1934, 
1935),  and  Deraniyagala  (1939).  More  recent  studies  on  this 


Taylor  (] 


vSsquez  Ruesta  (1982-83),  Webb  et  al.  (1983b),  Ferguson 
(1985),  Hall  (1985b),  and  Deeming  and  Ferguson  (1990). 


proportional  gross  morphology  of  extant  crocodilians  as  did 


1990),  Crocodvlus  niloticus  (Cott,  1961;  Graham,  1968,  1976; 


Crocodvlus  siamensis  (Muller,  1923,  1927),  and  Paleosuchus 
trigonatus  (Magnusson  and  Lima,  in  press) . The  works  of 
Parker  (1882) , Mook  (1921a  and  b) , Kalin  (1933) , Iordansky 
(1973),  Langston  (1973),  and  Greer  (1974)  were  confined  to 
cranial  growth  in  crocodilians.  However,  the  quantitative 


1955) , Brongersma  (1941) , and  Iordansky  (1973) . However, 
these  authors  made  no  quantitative  assessment  of  growth  with 
size,  sex,  or  age.  Dodson  (1975)  later  presented  allometric 
coefficients  for  27  attributes  of  cranial  growth  for  &. 
mississippjensis.  irrespective  of  sex  or  age.  Subsequently, 
Webb  et  al.  (1978)  reported  the  interorbital  width  to  be  a 
sexually  dimorphic  character  in  Australian  £.  oorosus . 
However,  the  predictive  values  of  their  discriminant  models 
to  correctly  identify  gender  were  low  except  for  larger 
animals;  an  evident  reflection  of  size  dimorphism. 

Similarly,  Montague  (1984b)  reported  sexual  dimorphism  in 
the  midpoint  width  of  the  cranial  roof  in  £.  novaeouineae. 
but  concluded  that  this  trait  was  an  unreliable  indicator  of 

study  Area  and  Methods 

Morphometry 

Morphometric  measurements  of  425  £.  novaeouineae  (342 
skulls  and  83  live  specimens)  were  recorded  during 
1980-1982.  These  include  162  males  and  83  females  plus  180 
individuals  of  unknown  gender.  All  animals  were  referable 
to  the  southern  form  of  this  species  (Hall,  1989) . Most 
specimens  were  collected  in  Lake  Murray  District,  Western 
Province,  PNG  (Fig.  5-3),  while  a handful  came  from  Central 
Province  (Brown  River,  Waigani  Swamp) , PNG.  Lake  Murray 
District  is  primarily  tropical  moist  forest  with  pronounced 
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e Fly  (left)  and  L 


wet  (Novenber-May)  and  dry  (June-Ootober)  seasons.  Annual 
mean  precipitation  is  306  cm  and  325  cm  in  Baboa  and  Pangoa, 
respectively,  in  the  Lake  Murray  subdistrict  and  258  cm  at 
Bosset  in  the  Middle  Fly  subdistrict  (Hall,  CHAPTER  3). 
Detailed  descriptions  of  the  study  area  are  presented 
elsewhere  (Hall,  1983;  Hall  and  Johnson,  1987). 

1 obtained  data  from  skulls  provided  by  hunters  and 
from  culled  captive  stock  caught  from  the  wild. 

Additionally,  I examined  two  skulls  from  Western  Province  at 
the  British  Museum  (Natural  History) , London.  All  skulls 
were  dried  at  ambient  air  temperature  for  at  least  30  days 
prior  to  inspection.  Pre-  and  post-drying  measurements  were 
made  on  culled  captive  animals  to  estimate  shrinkage  of 
certain  characters.  Thus,  correction  factors  can  be  applied 
to  data  that  I also  obtained  from  live  individuals  or  from 
preserved,  but  uncleaned  skulls. 

The  morphometric  measurements  taken  (Figs.  5-4  - 5-8), 
their  abbreviations,  explanation,  and  the  nearest  whole  unit 
of  value  employed  are  listed  in  Table  5-1.  Measurements  up 
to  150  mm  were  taken  with  0.05  mm  dial  calipers.  A steel 
tape  was  used  for  larger  dimensions.  Body  mass  was  weighed 
on  a Salter  clock  scale.  Manual  probing  of  the  cloaca 
(Chabreck,  1963)  was  used  to  determine  sex. 

Ratio  indices  were  used  for  skull  characters  to  reflect 
relative  growth.  The  indices  used  and  their  methods  of 
calculation  are  indicated  in  Table  5-2.  Cranial  indices 
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used  in  this  study  were  modified  from  Iordansky  (1973)  by 
substituting  dorsal  cranial  length  (DCL)  for  cranial  length 
(CL)  as  an  index  divisor.  The  latter  measurement  includes 
the  occipital  condyle  which  is  often  broken  during 
collection.  Also,  the  former  measurement  is  preferable 
since  it  can  be  recorded  from  living  individuals,  whereas 
the  latter  cannot.  Mandibular  indices  were  of  my  own 

Voucher  specimens  that  were  imported  to  the  United 
States  were  exported  from  PNG  under  the  following  CITES 
permit  numbers:  810057-810059,  and  810296-810299. 

Analyses 

Linear  regression  analysis  was  used  to  derive  equations 
for  the  estimation  of  DCL  from  body  size  measurements  or 
other  skull  attributes  used  as  independent  variables. 

Slopes  and  intercepts  of  regression  formulae  were  compared 
by  F-tests  to  evaluate  for  gender  specificity.  Coefficients 
were  also  calculated  for  the  prediction  of  body  size 
attributes  from  DCL  by  the  same  method. 

and  females  to  age  13  (Hall,  CHAPTER  3)  were  fitted  to  a 
reparameterized  Richards  curve  (Brisbin  et  al.,  1986)  of  the 
form:  w,  = [A<'-rt-(A<,'w-W0,'-",,)exp(-2t/T(m+l) ) ]'''" 
where  the  parameters  are  as  follows.  Wt  and  W„  represent 
known  body  mass  or  length  at  time  t and  at  time  t=0 


(hatching) ; A is  the  final  asymptotic  length  or  mass,  while 
T and  m represent  the  length  of  the  total  growing  period  and 
the  Richards  curve  shape  parameter,  respectively.  In  this 
study  Wt  is  335  cm  for  males  and  265  cm  for  females  (below) , 
while  A is  28  cm  for  both  sexes  (Hall,  1985b) . A second 
pair  of  models  was  constructed  with  respective  Wt  values  of 
360  and  285  (see  RESULTS) . The  parameter  T was  estimated  at 
50  years  for  males  and  females  and  m was  set  at  a value 
0.667  after  Brisbin  (1990).  A third  pair  of  models  using 
another  Richards  derivation,  the  Brody  curve  (Gates  and 
Dahm,  1989)  where  the  shape  parameter  m is  0,  was  calculated 
for  comparative  purposes  using  the  first  set  of  Wt  values. 
This  curve  takes  the  form:  W,  = A[l-Bexpe(w>] , where  W,  and  A 
are  as  above,  and  B and  C are  constants  with  values  of  0.9 
and  -0.2  in  the  models  presented. 

Gender  analysis  of  relative  growth  ratios  for  £. 
novaeouineae  was  tested  by  two  methodologies  to  assess  the 
comparative  robustness  of  each  approach  and  their  utility 
for  properly  classifying  observations.  A traditional 
parametric  multivariate  method,  discriminant  analysis 
(DISCRIM;  SAS,  1988) , was  first  used  followed  by  a 
non-parametric  binary  classification  tree  analysis  program 
(CART  [California  Statistical  Software,  Inc.];  Breiman  et 
al.,  1984).  The  CART  methodology  offers  several  important 
advantages  over  discriminant  analysis  in  that  it  treats 


observations  with  missing  values  and  that  it  states  the 
decision  rules  by  which  classification  trees  are  made. 

Both  DXSCRIM  and  CART  evaluate  classification  criteria 
by  their  performance  in  the  classification  of  future 
observations.  Initial  error  count  estimates  for  each  method 
were  derived  by  summing  the  number  of  misclassified  sex 
observations  from  the  test  data  sets  and  dividing  by  the 
total  number  of  observations.  The  resulting  error  rates  are 
optimistically  biased  and  referred  to  as  'apparent1  error 
rates  (AERs) . To  evaluate  the  correctness  of  these 
classification  rates,  crossvalidation  data  sets  were  created 
then  analyzed  for  each  method  to  reduce  model  biases  and 
associated  variances  and  to  establish  'true'  error  rates 
(TERs) . Crossvalidation  is  a form  of  jacknifing  that  treats 
n-l  out  of  n training  observations  as  a learning  set;  it 
derives  discriminant  functions  based  on  these  n-l 


observations  and  applies  them  to  classify  the  one 
observation  left  out.  Iterations  are  then  performed 
successively  for  each  of  the  remaining  n training 
observations  and  a misclassification  rate  is  calculated  as 


I further  analyzed  DISCRIM  predictor  models  by  using 
stepwise  discriminant  analysis  (STEPDISC;  SAS,  1988) . This 
procedure  selects  a subset  of  variables  from  the  original 
data  set  to  provide  a 'best  fit'  model  for  discrimination. 

, I then  used  the  character  variables  selected  by 


Similarly, 


CART  in  the  initial  analysis  above  to  create  a new  series  of 
comparative  'best  fit'  models  for  gender  classification. 


General  Morphology 

Linear  regression  coefficients  of  body  size  and  skull 
attributes  used  as  independent  variables  to  predict  dorsal 
cranial  length  (DCL)  are  listed  in  Table  5-3.  Most 
variables  showed  a slight  allometry  with  increasing  dorsal 
cranial  length  (Fig.  5-10) . Allometry  was  most  pronounced 
with  respect  to  body  mass  (BM;  Fig.  5-10d)  for  which 
logarithmic  transformation  (log10)  was  made.  Several 
cranial  (snout  length  [SL  - Fig.  5-10g] , cranial  length  [CL 
- Fig.  5-10e],  upper  ramus  length  [UL  - Fig.  5-10q])  and 
mandibular  measurements  (mandible  length  [ML  - Fig.  5-10v] , 
dentary  length  [DL  - Fig.  5-10y] , lower  ramus  length  [LM  - 
Fig.  5-10ab])  exhibited  near  perfect  (rz  > 0.99)  isometry  as 
predictors  of  DCL.  Correlations  of  body  size  and  skull 
measurements  with  DCL  were  uniformly  high,  often  extremely 
so.  The  fit  of  regression  coefficients  ( r2  values)  was 
consistently  better  for  males  than  females.  Variability  was 
most  pronounced  with  respect  to  interorbital  width  (IOW)  and 


dental  alveoli. 
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Predicting  Body  Size  From  DCL 


body  size  measurements  may  be  accomplished  by  the  regression 
coefficients  in  Table  5-4.  The  observed  shrinkage  between 
freshly  killed  C . novaeauineae  and  their  cleaned,  dried 
skulls  was  approximately  4%  DCL,  which  conforms  with  earlier 
reports  (Montague,  1984b) . Thus,  a correction  factor  of 
1.04  DCL  should  be  used  for  dried  skulls  prior  to  estimating 

If  DCL  is  an  estimated  value  derived  from  one  of  the 
preceding  isometric  predictors  (SL,  CL,  UM,  ML,  DL,  LM) , the 
standard  error  of  the  final  regression  estimate  (Ej)  is 
calculated  as  follows:  Jt  (b2E,)J+E28) , where  b2  is  the  slope 
of  the  second  equation,  and  E1  and  E2  are  the  error 
estimates  from  the  first  and  second  equations,  respectively 
(Webb  and  Messel,  1978). 

The  largest  live  male  and  female  in  this  study  measured 
313  and  265  cm  total  length  (TTL) , respectively:  while  the 
largest  male  skull  measured  455  mm  DCL  (473  mm  pre-shrinkage 
estimate)  which  would  equate  to  a TTL  estimate  of  360  cm. 

The  maximum  verified  sizes  (TTL)  of  southern  form  C. 
novaeauineae  are  335  cm  for  males  (Montague,  1984b)  and  265 
cm  for  females  (above) . The  form  of  this  species  that 
occurs  in  northern  New  Guinea  attains  larger  known  TTLs: 
males  - 347  cm;  females  = 305  cm  (Hall,  CHAPTER  3). 
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Sexual  Dimorphism 


Sexual  dimorphism  was  observed  in  snout-vent  length 
(SVL) , TTL,  and  in  21  o£  32  skull  attributes  (Table  5-5). 
These  differences  varied  from  slight  to  pronounced  at  the 
size  of  the  largest  female  in  this  study  (352  mm  DCL) . Male 

cranial  length  (CL) , and  mandible  length  (ML) , but  have  a 
shorter  premaxillary  symphysis  length  (PXS)  and  upper  ramus 
length  (UM)  than  females.  Males  also  exhibit  greater  widths 
with  respect  to  their  cranium  (CW) , snout  (SH) , orbit  (OH) , 
interorbit  distance  (IOH) , cranial  roof  (HCR) , external 
nares  (WN) , mandibular  symphysis  (WMS) , and  distance  across 
the  surangulars  (HSR) . However,  the  most  obvious  evidence 
of  dimorphism  exists  in  regard  to  dentition  in  which  the 
alveoli  PM4,  M5,  Dl,  D4,  and  Dll  are  conspicuously  more 
robust  in  males.  The  alveoli  PM4,  M5,  Dl,  and  D4  are 
pseudo-canines  in  the  genus  Crocodvlus.  Their  dimorphic 

the  social  ordering  of  adult  males  that  compete  for 
territories  and  access  to  breeding  females. 

Changes  in  the  directional  rates  of  skull  growth 
indicative  of  sexual  dimorphism  commence  in  Q.  novaeauineae 
at  about  100  cm  TTL.  Females  begin  to  exhibit 
proportionally  longer  skulls  than  males  of  equivalent  TTL; 
however,  differences  in  the  proportional  widths  of  skulls 
are  not  manifest  until  a TTL  of  about  180  cm  is  attained. 


Isslpl8"sl|s^||ssas3s 


iig,;^.;gnjj^g.g:  s^'s^rk  «,. 


Il:i 


!S:i 


S«:i 

’ll 

S:S 


3:J 


*S:| 

S:1 


!:S" 

S:SS5 

!:” 

!:S 


Size-Age  Relations 


Growth  estimates  for  male  £.  novaeouineae  to  age  nine 
and  females  to  age  13  (Hall,  CHAPTER  3)  were  fitted  to  a 
modified  Richards  curve  (Brisbin  et  al.,  1986)  to  derive 
growth  models  for  predicting  size-age  relations.  Initial 
models  were  calculated  from  known  maximum  sizes  of  335  and 
265  cm  TTL  for  males  and  females,  respectively,  as  final 
asymptote  size  parameters.  These  equations  were  rerun  to 
develop  a second  model  set  using  the  predicted  maximum  size 
estimate  of  360  cm  TTL  for  males  (above)  and  with  a value  of 
285  cm  TTL  for  females.  This  latter  value  was  obtained  by 
using  a proportional  size  ratio  estimate  as  for  known  sizes 
(265/335  x 360) . A third  set  of  growth  models  were 
calculated  using  a Brody  curve  (Gates  and  Dahm,  1989) , 
another  modification  of  Richards  curve,  with  final  asymptote 
sizes  as  in  the  first  model  pair  (Fig.  5-11) . 

Although  the  Brody  growth  model  closely  paralleled  the 
female  model  with  known  asymptote  size,  it  considerably 
underestimated  the  comparable  model  for  males.  Regardless 
of  gender,  growth  forms  diverged  at  about  age  15  for  all 
models.  Although  these  models  may  approximate  mean  growth 
for  wild  £.  novaeauineae.  they  are  speculative  and  their 
accuracy  is  unknown.  However,  it  is  apparent  that 
considerable  error  could  result  from  attempts  to  age  animals 
longer  than  210  cm  that  are  beyond  the  point  of  model 
divergences. 


The  original  relative  growth  data  set  incorporated  a 
total  of  27  ratio  indices  for  a complete  skull.  These 
included  15  cranial  and  12  mandibular  aspects  of  growth 


There  was  considerable  overlap  among  canonical 
coefficients  in  the  relative  growth  data  set.  Nevertheless, 
initial  discriminant  analysis  (DISCRIM)  results  were  highly 
promising  (e.g.,  overall  apparent  error  rate  [AER]  = 3.6% 
misclassification  with  a complete  skull) , but  alarmingly 
illusory  as  the  crossvalidation  true  error  rate  (TER)  was 


44.8%.  When  the  same  data  set  was  evaluated  by 
classification  tree  analysis  (CART)  the  resulting  AER 
(26.1%)  provided  a realistic  estimate  of  the  TER  (28.7%). 

The  best  overall  discrimination  resulted  from  CART  analysis 
of  the  cranium  only,  where  a subset  of  only  two  cranial 
attributes  (relative  width  of  snout  [RWST]  and  relative 
length  of  cranial  roof  [RLR])  were  selected  by  the  model  for 
classification  of  gender  (Fig.  5-12) . 

The  wide  disparity  between  AERs  and  TERs  was  a 
consistent  feature  of  DISCRIM,  while  CART  results  yielded 
error  rates  that  much  more  closely  paralleled  each  other. 
Whether  a complete  or  partial  skull  was  analyzed,  fewer 
observations  were  misclassified  using  CART  methodology  than 
when  TERs  are  taken  into  account.  Regardless 
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of  the  method  used,  better  gender  classification  resulted 
from  use  of  the  cranium  than  the  mandible  alone. 

subsequent  examination  of  the  original  DISCRIM  data  set 
by  stepwise  discrimination  analysis  (STEPDISC)  resulted  in 
reduced  model  sizes  consisting  of  a maximum  of  four 
variables  (RWST,  RLE,  relative  upper  ramus  length  [RLUMR], 
relative  mandibular  symphysis  width  [ RWSS ] (Table  5-7) . 
Construction  of  a series  of  'best  fit'  models  by  DISCRIM 
using  the  STEPDISC  selected  variables  gave  much  improved 
results  and  better  correlation  among  AERs  and  TERs  that 
approximated  the  CART  'best  fit'  model  consisting  of  the 
variables  RWST  and  RLR  (Table  5-8) . 

Relative  growth  Patterns 

Relative  growth  patterns  for  the  skull  of  Crocodvlus 
novaeauineae  are  shown  in  Figure  5-13.  Growth  is 
characterized  by  three  general  stages:  (I)  a hatchling 
phase,  (II)  a juvenile  through  middle-aged  adult  phase,  and 
(III)  an  old-aged  adult  phase. 

Stage  I is  maintained  through  most  of  the  first  year  of 
growth  until  a DCL  of  60-70  mm  (40-50  cm  TTL)  is  attained. 

In  this  stage,  the  preorbital  region  of  the  skull  is 
characterized  by  a relatively  short,  robust  snout  with 
comparatively  narrow  external  nares.  The  orbital  region 
features  enormously  enlarged  orbits  (primarily  in  relative 
width)  with  correspondingly  narrow  interorbital  width;  while 
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in  the  postorbital  portion  of  the  skull,  the  cranial  roof 

The  change  from  stage  1 to  stage  II  skull  growth  is 
very  abrupt,  paralleling  the  'saltatory  ontogeny'  model 
proposed  by  Balon  (1984) . The  secondary  stage  is  a 
protracted  phase  exhibited  until  animals  reach  a DCL  of 
about  350  mm  (roughly  260  cm  TTL) . It  is  in  this  stage  t 
the  crocodilian  skull  acquires  the  proportions  of 
generalized  morphology  typically  associated  with  the 
species.  The  snout  becomes  relatively  long  and  slender, 
*idth  in  external  nares  and 


interorbital  distance,  the  orbits  are  much  reduced  in  size 
from  stage  I,  and  the  cranial  roof  becomes  a concave-shaped 
trapezoid  for  the  duration  of  life. 

By  virtue  of  size,  stage  III  growth  consists 
predominantly  of  large  adult  males  and  excludes  all  but  the 
largest  size  females.  In  stage  III  the  directional  rates  of 
growth  undergo  another  transformation.  As  the  skull  length 
asymptote  is  approached,  a reversal  toward  the  hatchling 
stage  I phase  is  manifest  with  respect  to  the  skull  becoming 
increasingly  more  robust  while  the  snout  becomes 
proportionally  shorter.  Thus,  an  individual  may  exhibit 
identical  relative  growth  ratios  for  a given  attribute 
within  its  life  span,  but  at  different  stages  of  growth  and 
at  widely  differing  ages.  Hence,  some  independent  measure 


of  size,  such  as  DCL,  must  be  incorporated  to  lend 
meaningful  interpretation  to  relative  growth  data. 


Relative  crewth 

Relative  growth  indices  are  a useful  way  to  chronicle 
ontogenetic  patterns  of  growth  and  to  segregate 
morphologically  similar  species  (Dodson,  1975;  Simpson, 

1980;  Hall,  1989).  Such  indices  in  this  study  reveal  three 
distinct  stages  of  skull  growth  in  £.  novaegulneae:  a 
hatchling  phase  (I) , a juvenile  to  middle-aged  adult  phase 
(11) , and  an  old-aged  adult  phase  (III) . 

Saltatory  ontogeny  (sensu  Salon,  1984)  is  a series  of 
near  spontaneous  changes  in  form  and  function  that  lead  to 
'threshold'  states  between  the  organism  and  its  environment. 
This  life-history  model  is  a hierarchical  system  of 
intervals  based  on  evolutionary  opposite  trends  between 
altricial  (generalist)  and  precocial  (specialist)  forms. 
Examples  include  threshold  responses  from  photophobia  to 
phototaxia  and  from  endogenous  to  exogenous  feeding  in  early 
ontogeny.  The  abrupt  shift  in  patterns  exhibited  by 
hatchling  £.  novaeguineae  from  stage  I to  stage  II  relative 
skull  growth  mimics  a saltatory  response.  Hatchling 
crocodilians  possess  a residual  embryonic  yolk  reserve  that 
continues  to  provide  nutrients  for  growth  in  the  early 


posthatchling  period  (Whitehead,  1990) . 


addition 


depletion  of  yolk  reserves  during  the  stage  I phase  of 
hatchling  cranial  growth,  there  is  a concomitant  dietary 
shift  from  an  invertebrate  prey  base  toward  a gradual 
inclusion  of  vertebrate  items  that  form  the  mainstay  of 
forage  for  older  animals.  This  is  evidenced  by  structural 
modifications  of  the  skull  whereby  hatchlings  are 
principally  visually  oriented  versus  the  diminution  of  the 
orbital  region  and  expansion  of  the  external  nares 
throughout  growth  stages  II  and  III,  inferring  greater 
reliance  on  olfactory  prowess  in  older  animals  rather  than 
on  visual  acuity.  The  shift  from  stage  II  to  stage  III 
skull  growth  is  not  as  rapid  and  pronounced  as  from  stage  I 
to  stage  II,  but  it  is  in  this  stage  that  animals  attain  the 
size  and  are  afforded  the  massiveness  and  strengthening  of 
the  skull  through  increased  robustness  to  feed  on  the 
largest  vertebrate  prey  items  available. 


gender  Classification 

The  correctness  of  classification  rates  of  biological 
observations  by  discriminant  analysis  are  strongly 
influenced  by  group  sample  sizes  and  the  distributional 
overlap  among  groups.  When  overlap  is  high,  as  in  this 
study,  sample  sizes  must  be  increased  with  respect  to  the 
number  of  variables  to  provide  a given  level  of  precision. 
Williams  et  al.  (1990)  recommend  that  sample  sizes  be  at 


variables 


measured 


among  group  overlap  is  low  and  substantially  more  as  the 
overlap  increases.  Sample  sizes  in  this  study  for  a 
complete  skull  ranged  from  approximately  four  times  the 
variable  number  for  the  general  DISCRIM  model  (Table  5-6)  to 

As  noted  earlier,  CART  methodology  maximizes  the  number  of 
available  observations  since  it  treats  missing  value  data, 
whereas  DISCRIM  does  not.  This  feature,  plus  its  ease  of 
use  and  interpretation,  makes  CART  preferable  to  DISCRIM  for 
classification  purposes.  More  importantly,  in  this  study 
CART  provided  consistently  lower  TERs  that  gave  more 
realistic  assessments  of  model  AERs  than  did  DISCRIM. 
However,  the  choice  of  methodology  is  not  so  clear  cut  when 
one  compares  'best  fit'  models  since  DISCRIM  performed 
nearly  as  well  as  CART  in  that  regard,  although  it  entailed 
the  use  of  three  more  attributes  (IMS,  UM,  WMS)  than  the 
latter  method  (both  used  DCL,  LCR,  SL,  SW  and  WCR) . 

The  fact  that  the  cranial  roof  and  snout  were  selected 
as  predictive  variables  for  the  classification  of  gender  in 
both  the  DISCRIM  and  CART  'best  fit'  models  suggests  that 
better  results  might  be  obtained  through  the  addition  of 
other  measures  of  the  shapes  of  those  variables.  Post-  and 
preorbital  attributes  that  may  prove  of  use  in  further 
investigation  are  as  follows:  (1)  the  midpoint  width  of  the 
cranial  roof,  and  (2)  the  intersuprafenestral  width  of  the 
parietal:  plus,  three  other  measurements  of  the  snout. 


These  latter  include  (3)  the  mid-snout  width  (at  the  level 
of  the  5th  maxillary  alveoli),  (4)  the  width  at  the  point  of 
anterior  snout  constriction  (posterior  to  the  external 
nares) , and  (5)  the  maximum  premaxillary  snout  width  (at  the 
midpoint  of  the  external  nares) . 


Sise-Ace  Relations 

Attempts  at  aging  crocodilians  usually  involve 
capture-recapture  studies  (Mcllhenny,  1934,  1935;  Lawson, 
1972;  Staton  and  Dixon,  1975;  Webb,  1977;  Gorzula,  1978; 

Webb  and  Messel,  1978;  Webb  et  al.,  1978;  Chabreck  and 
Joanen,  1979)  or  observed  growth  in  captivity  (Cott,  1961, 
Montague,  1982) . Graham  (1968)  employed  growth  annuli  in 
the  dentary  bone  and  the  weight  of  the  eye  lens  to  determine 
the  age  of  a sample  of  harvested  Q.  niloticus.  However, 
those  methods  provided  less  than  satisfactory  results. 
Subsequently,  Hutton  (1986,  1987c)  reported  that 

better  resolution  and  a more  reliable  method  of  aging. 
However,  as  noted  above,  considerable  margin  for  error  in 
aging  crocodilians  still  exists.  Hutton's  (1986)  study 
provides  the  most  rigorous  evidence  to  date  of  the  magnitude 
of  error  resulting  from  attempts  to  age  these  reptiles.  In 
that  study,  the  range  of  errors  of  individual  e 

errors  of  mean  estimates  ranged  from  15%  f 


year  old 


animals  to  9t  for  a 46  year  old  individual.  Errors  are 
likely  to  be  greater  in  females  due  to  lines  of  arrested 
growth  resulting  from  reabsorption  of  calcium  from  ossified 
areas  during  the  reproductive  cycle  (Wink  and  Elsey,  1986) . 
Thus,  use  of  the  growth  curves  in  this  study  (Fig.  5-11) 
will  likely  produce  spurious  results  for  all  but  the 
youngest  age  classes. 

The  teeth  of  certain  mammalian  species  exhibit  growth 
annuli  and  have  long  been  used  for  aging  purposes  (e.g.. 
Marsh,  1981) . Crocodilian  teeth  also  possess  growth  zones 
(Poole,  1961;  Johnston,  1979) . However,  crocodilians  have 
polyphyodont  dentition  characterized  by  replacement  in 
irregular  seguence  and  at  irregular  intervals  (Edmund, 

1962) . These  properties  render  crocodilian  teeth  unsuitable 
for  aging  purposes.  I believe,  though,  that  the  size  of 
dental  alveoli  may  prove  to  be  a useful  indicator  of  age  in 
crocodilians  since  growth  of  an  alveolus  appears  to  be 
continual,  except  in  rare  pathological  conditions  (Hall, 
1985a) . Further  investigation  along  these  lines  is 
required. 

Summary 

derived  with  a high  degree  of  accuracy  from  a complete  skull 
or  from  fragments  thereof.  Allometric  coefficients  are 
provided  for  the  estimation  of  dorsal  cranial  length  (DCL) 


from  body  and  skull  attributes  and  for  the  determination  of 
body  size  attributes  from  DCL.  Best  estimates  of  DCL  are 
obtained  from  the  use  of  cranial  attributes  SL  (snout 
length) , CL  (cranial  length) , UM  (upper  ramus  length)  and 
mandibular  attributes  ML  (mandible  length) , DL  (dentary 
length),  LM  (lower  ramus  length). 

Sexual  dimorphism  is  subtle  but  significant  with 
respect  to  body  size  (SVL,  TTL)  and  21  of  32  skull 
attributes.  Differences  among  sexes  ranged  from  1.3%  to 
28.3%  for  a given  attribute  at  the  size  of  the  largest  known 
female  £.  novaequinsas  and  were  most  pronounced  in  the 
relative  sizes  of  dental  alveoli. 

Relative  growth  indices  reveal  three  distinct  stages  of 
skull  growth.  Stage  I is  characteristic  of  hatchling 
crocodiles,  stage  II  of  juvenile  through  middle-aged  adults, 
and  stage  III  of  large  elderly  adults  (principally  males) . 
Stage  I growth  abruptly  shifts  to  stage  II  at  about  60-70  mm 
DCL  and  mimics  a saltatory  ontogenetic  response.  The  shift 
from  stage  II  to  stage  III  occurs  at  about  350  mm  DCL  and  is 
less  abrupt  than  stage  I.  Stage  III  is  characterized  by  a 
reversal  in  directional  rate  of  growth  toward  that  found  in 
stage  I where  the  relative  snout  length  is  decreased  and 


The  use  of  relative  growth  indices  for  the  skull  of  £. 
novaeouineae  provides  a reasonable  means  of  determining  the 
gender  of  unknown  sex  animals  by  classification  through 


multivariate  discriminant  analysis  (DISCRIM)  or  through 


non-parametric  binary  tree  classification  (CART) . CART 
methodology  provided  consistently  lower  apparent  error  rates 
(AERs)  that  better  correlated  with  crossvalidation  true 
error  rates  (TERs)  than  those  obtained  by  DISCRIM.  However, 
when  'best  fit'  models  were  selected  from  variables  used  in 
general  classification  models,  DISCRIM  performed  nearly  as 
well  as  CART.  Error  rates  were  invariably  greater  for 
females  than  males  regardless  of  the  method  or  model 
employed.  Suggestions  are  made  for  the  possible  improvement 
of  classification  results  through  the  addition  of  other 
measures  of  predictive  attributes  used  by  the  'best  fit' 


conclusion,  the  use  of  hunter  harvested  skulls  of  £. 
provide  valuable  demographic  insights  to 
ix  composition  of  populations  from 
locations  that  are  otherwise 
This  knowledge  will  offer  a 
monitoring  population  trends  of 


remote,  difficult  access 
unlikely  to  be  obtained, 
cost-effective  method  01 
future  harvests  and  enh; 
achieving  optimal 


CHAPTER  6 
CONCLUSION 


The  potential  for  drawing  demographic  inferences  of 
crocodilian  populations  from  seemingly  disparate  sources  and 
on  hitherto  poorly  known  aspects  of  their  biology  was 
investigated  in  this  study.  The  patterns  of  harvest, 
reproduction,  and  skull  morphometry  were  examined  variously 
for  one  or  more  of  three  species:  populations  of  New  Guinea 
(Crocodvlus  novaequineae)  and  saltwater  (£.  oorosus) 
crocodiles  from  Papua  New  Guinea  (PNG) , plus  populations  of 
American  alligators  ( Alligator  raississipoiensisl  from  the 
states  of  Florida  and  Louisiana,  USA. 

Historically,  tens  of  millions  of  crocodilians  have 
been  killed  as  vermin  or  for  their  commercial  value.  Annual 
worldwide  trade  is  presently  estimated  to  be  between  one  and 
two  million  hides,  of  which  perhaps  70%  or  more  are  illegal 
(Brazaitls,  1987).  Hence,  accurate  harvest  data  are  seldom 
available  for  interpretation  and  such  data  as  exist  often 
are  based  on  small  samples  that  result  from  the  interception 
of  confiscated  hides  (e.g.,  Rebelo  and  Magnusson,  1983).  A 
notable  exception  is  the  study  of  Taylor  and  Neal  (1984) 


Louisiana  alligator 


My  analysis  of  1969-1980  harvest  data  for  PNG  crocodile 
populations  (CHAPTERS  2 and  3)  indicated  that  substantial 
numbers  of  both  species  were  continuing  to  be  harvested  on 
an  apparent  sustained  yield  basis  in  spite  of  several 
decades  of  previously  unregulated  exploitation.  Although 
management  controls  were  instituted  in  1969,  the  crocodile 
skin  trade  remained  largely  in  the  hands  of  expatriates 
until  independence  was  attained  in  1975.  Following  that, 
collaborative  efforts  were  made  with  the  Food  and 
Agricultural  Organization  of  the  United  Nations  (FAO)  to 
restructure  the  industry  so  that  indigenous  landowners 
became  the  focal  beneficiaries  of  the  resource  (Bolton  and 
Laufa,  1982) . During  these  years  of  regulated  harvest  a 
tremendously  large  harvest  data  base  was  developed; 
unfortunately,  the  compilation  was  lacking  and  much  was  left 
to  desire  in  the  type  of  data  generated.  Nonetheless,  it 
was  possible  for  me  to  reconstruct  species-specific  harvest 
size  and  sex  compositions  on  those  data  to  gain  demographic 
insights  to  the  structures  of  previous  populations.  While 
this  entailed  the  use  of  some  rather  liberal  assumptions, 
the  process  also  afforded  the  opportunity  to  rectify  data 
collection  procedures  for  future  harvests  to  provide 
management  authorities  with  a more  instructive  format  on 
which  to  base  regulatory  decisions.  Those  reforms  were 
instituted  and  subsequent  monitoring  indicates  that 


population  recovery  of  these  species  is  continuing 
satisfactorily  (Hollands,  1987) . 

Even  for  the  best  studied  crocodilian  species,  A • 
mississipniensis.  the  understanding  of  demographic 
parameters  often  amounts  to  blatant  speculation 
(Abercrombie,  1989) . Speculative  or  not,  the  analysis  in 
CHAPTER  3 offers  a method  of  estimating  cumulative  and 
age-specific  mortality  rates  from  harvest  returns  and 
provides  some  insight  as  to  what  may  be  occurring  (e.g., 
greatly  higher  mortality  among  harvested  populations  of  C. 
porosus  than  £.  novaeauineael . Similarly,  CHAPTER  4 
outlines  a simple  but  effective  manner  of  accurately 
determining  structures  of  nesting  female  crocodilian 
population  segments.  The  inferential  evidence  on  the 
effects  of  density-dependent  nesting  as  a consequence  of 
adult  female  size  offers  an  intriguing  hypothesis  for  the 
type  of  potential  experimental  manipulation  advocated  by 
Macnab  (1983,  1985)  for  wildlife  managers. 

However,  there  is  a more  immediately  applicable  lesson 
from  the  demographic  inferences  of  CHAPTER  4 . This  is  that 
there  is  yet  another  line  of  independent  evidence  (see 
Montague,  1984 1 Cox,  1985,-  Hall  and  Johnson,  1987) 

brood  stocks  in  PNG  is  less  effective  than  originally 
intended.  That  the  original  act  was  drafted  without  benefit 
of  adequate  biological  input  is  an  important  consideration. 


In  the  mid-1980's  a crocodile  management  program  modeled 
after  the  one  designed  in  PNG  was  established  in  neighboring 
Irian  Jaya  (Cox,  1990) . Having  benefitted  from  the  PNG 
experience,  slot  size  limits  on  the  trading  of  Irian  Jayan 
crocodile  skins  were  set  at  25  to  46  cm  commercial  belly 
width  (BW)  as  opposed  to  limits  of  18  to  51  cm  in  PNG. 
Although  not  wholly  effective,  a 46  cm  BW  size  limit  (101  cm 
snout-vent  length  (SVL)  equivalent)  affords  greater 
protection  to  nesting  female  £.  novaeauineae  than  does  a 51 
cm  size  limit  (see  Fig.  4-3) . Serious  consideration  should 
be  given  to  revising  the  slot  size  limits  in  PNG  and  to 
standardizing  them  with  those  in  Irian  Jaya.  This  would 
help  thwart  the  possible  transhipment  of  undersized  or 
oversized  skins  from  Irian  Jaya  to  PNG  and  simultaneously 
better  fulfill  the  goals  of  the  initial  legislation  on  upper 
size  trading  limits. 

The  options  available  to  management  authorities 
responsible  for  the  conservation  of  these  resources  may  be 
strikingly  different  for  as  much  economic  and  logistic 
reasons  as  for  biological  ones.  The  methodologies  I offer 
for  the  assessment  of  demographic  parameters  of  crocodilian 
populations  should  prove  cost-effective  for  developed  and 
developing  nations  alike.  The  use  of  hunter  harvested  (or 
illegal  kill)  skulls  as  a valuable  demographic  artifact 
(CHAPTER  5)  will  probably  prove  of  greater  benefit  to 
developing  tropical  nations  that  are  more  constrained  by 


logistic  and  economic  considerations.  An  effective  and 
inexpensive  sampling  program  could  easily  be  established  in 
remote  areas  subject  to  hunter  harvest  for  the  collection  of 
skull  parts  from  culled  animals.  A skull  provides  a much 
better  estimate  of  live  animal  size  (SVL  or  total  length) 


a minimum,  a single 


jaw,  the  dentary. 


of  harvests.  A sex  ratio  interpretation  for  C.  novaeguineae 
harvests  may  also  be  inferred  with  a reasonable  degree  of 
error,  providing  that  the  cranium  is  used  for  gender 
classification.  Further  study  on  this  subject  is  required 
to  refine  the  technique  and  to  assess  its  applicability 
toward  other  species. 

Following  the  termination  of  FAO  input  to  the  PNG 
crocodile  skin  industry  much  financial  support  for  the 


:he  program  management  has  been  derived  from 
private  sector  commercial  beneficiaries  to  supplement 

from  government  sources.  Whether  there  or  in  the 
s appropriate  and  vital  that  such  contributions  be 
form  of  'value-added'  conservation  (Hines  and 


Percival,  1987).  It  is  equally  imperative  that  management 
agencies  do  not  become  fiscally  undermined  by  governmental 
budgetary  processes  or  compromised  as  a consequence  of  such 
sources  and  that  they  be  able  to  maintain  active  and 


continual  research  components 


conservation  planning. 
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Philip  M.  Hall  was  born  on  13  June  1949  in  Jamaica,  New 
York,  and  subsequently  raised  in  Levittown,  New  York. 

Trying  not  to  let  school  interfere  with  his  education,  he 
left  Hofstra  University  after  one  year  in  the  elusive 
pursuit  of  some  grander  purpose.  Following  six  years  of 
Wander jahren  and  several  premature  mid-life  crises,  he 
returned  to  the  sphere  of  formal  education  at  North  Dakota 
State  University  in  1974.  He  received  an  A.A.S.  in  ecology 
and  wildlife  in  1975,  and  in  1976  graduated  with  a B.S.  in 
zoology. 

Shortly  thereafter  he  embarked  on  broadening  his 
cosmopolitan  perspectives  by  joining  the  U.S.  Peace  corps  as 
a volunteer  in  a Nepalese  fisheries  development  program. 
After  two  and  a half  years  of  introspective  irreverancy  and 
hitherto  unrealized  hardship  and  adventure  he  returned  to 
the  U.S. ; whereupon,  he  was  hopelessly  maladjusted  to  the 
ostentatious  ways  of  the  developed,  20th  century  Western 
World.  His  fears  of  remaining  a cultural  and  social  pariah 
were  soon  tempered  by  the  prospects  of  returning  to  the 
Third  (and  by  then,  normal)  World.  Following  enrollment  at 
the  University  of  Idaho  in  1979,  he  became  a United  Nations 
volunteer  and  went  to  Papua  New  Guinea  where  he  spent  two 


years  as  a crocodile  biologist  under  contract  to  the  Food 
and  Agricultural  organization  of  the  united  Nations. 
Afterward,  he  returned  to  complete  his  degree  program  and 
received  an  M.S.  in  zoology  from  Idaho  in  1983. 

In  the  aftermath  of  several  later  life  crises,  he 
continued  graduate  study  at  Louisiana  State  University 
during  1985-1986  while  intermittently  researching  alligators 
in  southwestern  Louisiana  and  working  aboard  Japanese 
trawlers  and  a Soviet  processor  in  the  Bering  Sea  and  Gulf 

In  January  1987  he  enrolled  at  the  University  of 
Florida  where  he  is  a candidate  for  the  degree  of  Doctor  of 
Philosophy,  in  the  Department  of  Wildlife  and  Range 
Sciences,  School  of  Forest  Resources  and  Conservation. 

During  his  years  he  has  had  the  good  fortune  to  have 
loved,  known,  and  worked  with  many  fine  people  and  is  fully 
appreciative  of  the  privileges  accorded  him  by  his  chance 
birthright  in  a democratically  free,  non 


i-totalitarian 


. rSP" 

t=— 


